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The  probability  that  part  of  a  set  of  equicorr elated  normal 
variables  are  positive  is  defined  by  a  multiple  integral  expression 
involving  the  multivariate  normal  density  function.   Although  much 
research  related  to  this  integral  expression  has  been  published,  most 
results  do  not  include  a  practical  method  of  its  evaluation.   Also, 
when  the  correlation  is  negative,  no  direct  method  of  evaluating  the 
integral  expression  is  available.   In  this  paper  we  discuss  several 
methods  of  expressing  the  integral.   One  of  these  expressions,  valid 
for  both  positive  and  negative  correlation,  is  nsed  to  obtain  numer- 
ical results. 

A  transformation  is  used  to  simplify  the  integral  expression 
for  the  probability  that  part  of  a  set  of  equicorr el a ted  normal  var- 
iable are  positive.   Then  the  probability  can  be  written  as  an  integral 
involving  the  real  normal  distribution  function  when  the  correlation 
is  positive,  and  the  complex  normal  distributions  function  when  the 
correlation  is  negative.   For  positive  correlation,  this  integral 
expression  has  been  used  by  other  authors  to  obtain  numerical  results. 

viii 


Next,  we  use  a  result  connecting  the  terms  of  a  binomial  series 
v/ith  Tchebycheff-Hermite  and  Legencire  polynomials  to  obtain  a  finite 
series  expression  for  the  probability.   Although  the  general  term  in 
the  series  involves  an  integral  which  cannot  be  evaluated  in  closed 
form,  this  integral  depeiids  only  on  the  correlation  and  can  be  eval- 
uated by  numerical  integration  for  both  positive  and  negative  correla- 
tion.  The  numerical  results  are  included  in  the  appendixes. 

As  an  application,  the  number  of  observations  larger  than  the 
sample  mean  is  used  as  a  test  statistic  for  testing  the  hypothesis 
that  a  population  is  normally  distributed.   The  small  sample  null 
distribution  is  derived  and  numerical  results  are  given.   Approxima- 
tions to  the  null  distribution  are  also  discussed. 

Finally,  we  discuss  three  other  methods  of  expressing  the 
probability  that  the  entire  set  of  equicorrelated  normal  variables  is 
positive.   Two  of  these  methods  express  the  probability  as  an  infinite 
series.   However,  in  both  cases  the  convergence  is  quite  slow.   The 
third  expression,  involving  moments  of  extreme  order  statistics,  can 
be  used  for  obtaining  numerical  results  only  for  limited  positive 
values  of  the  correlation. 


IX 


CHAPTER  1 
INTRODUCTION 

1. 1  Introduction 

In  a  recent  paper,  David  (1962)  suggested  using  the  number  of 
observations  larger  than  the  sample  mean  as  a  test  for  the  homogeneity 
of  a  random  sample.   Assuming  a  normal  population,  he  showed  that  the 
proportion  of  observations  larger  than  the  sample  mean  has  an  asymptotic 
normal  distribution.   However,  David  did  not  discuss  the  small  sample 
distribution  for  the  test  statistic. 

The  work  on  this  dissertation  began  in  search  for  the  small 
sample  distribution  for  David's  statistic.   However,  this  work  soon 
led  to  the  more  general  problem  of  finding  the  probability  that  part 
of  a  set  of  equicorrelated  normal  variables  are  positive  and,  in 
particular,  the  problem  of  evaluating  a  multivariate  normal  integral 
expression  for  the  probability  that  all  the  variables  are  positive. 

Much  research  related  to  the  multivariate  normal  integral  has 
been  published.   Gupta  (1963b) ,  in  addition  to  an  excellent  survey 
paper,  gives  a  complete  bibliography  of  articles  related  to  the  multi- 
variate normal  integral.   However,  only  a  few  of  these  articles  offer 
a  practical  method  of  evaluating  the  integral.   Also,  although  Steck 
(1962)  gives  a  relation  connecting  the  results  for  positive  and  nega- 
tive correlation,  no  direct  method  of  evaluating  the  integral  has 


been  obtained  when  the  correlation  is  negative.   Only  cwo  authors, 
Ruben  (1954)  and  Gupta  (1963a),  give  numerical  results. 

It  is  the  purpose  of  this  dissertation  to  find  at  least  one 
method  of  evaluating  the  multivariate  normal  integral  that  works  for 
both  positive  and  negative  correlation  and  that  can  be  used  easily  to 
obtain  numerical  results.   Then  the  small  sample  distribution  of  David's 
statistic  can  be  given  as  an  application  to  the  more  general  problem. 


1.2  Definition  of  P   (p ) 

r  :m 


Suppose  the  m  variates  X  ,X  , ... ,X  ,  each  with  zero  mean  and 

unit  variance,  have  a  multivariate  normal  distribution.   Of  interest 

is  the  probability  that  exactly  r  of  these  m  variables  are  positive. 

If  the  variables  are  mutually  independent,  the  problem  has  the 

binomial  solution   (  J  (—J  .   However,  when  the  variables  are  dependent, 

no  simple  solution  exists.   In  this  paper  we  shall  consider  the  case 

when  the  variables  have  common  correlation  o.  — r  <  0  <  1.   P    (p ) 

m-1  r:m 

will  denote  the  probability  that  exactly  r  of  the  m  variables  are 
positive.   That  is, 

P    (p)  =  EP(X±   >0 Xi  >  0;  X±    <0,...,Xi  <0), 

r:m  -"-1  xr       xr+l  m 

where  the  summation  is  over  all  partitions  fi  ,...,?.  ;  i   .  ....  ,i  1 

L  1      r   r+1      mJ 

of  the  set  fl,2,...,m}.   Since  X  ,X_,...,X  are  identically  distrib- 

J  12m 

uted,    the   above   equation  may   be  written    as 


P        (p  )    =    f"1  V  (X,   >  0 X     >  0;    X  <0 X     < 

r:m  ^  \r/  1  r  r+1  m 


Letting  g(x„ ,x  , .  .  .  ,x  )  represent  the  density  on  X„ ,X„, . . . ,X  , 

1   2      m  12m 

we  have 


0) 


d.2.1)         p     (P)  =  £)/.../  !■■■!  g(xi,x2>...xffi)dxi,dx2,...,dxm 

X  >  O  X  .  <  o 

i  i 

i<  r  i  >  r 


J   g(x1(x2(...,xm)  dxi(dx2(...>dxm, 


[r,m] 


where  [r,m]  will  be  used  to  denote  the  range  of  integration 

fx     >0:l<i<r;x     <0:    r+1  <  i  <  m]    . 
L    i  l  J 

In  the  case  r  =  m,  we  will,  for  convenience,  write 

P  (c  )  rather  than  P    (p  ) . 
m  H  m:m  r 

A  first  approach  to  the  problem  might  be  to  write  down  the 

density  function  g(x  ,x0 x  ).   From  multivariate  theory  we  have 

12m 

m  11      t-1 

(1.2.2)  gi\,*2 V    =    (2TT)         lVl         6 

_    oo   <    x       <    °° 

i 
i  =   1,2, . . . ,m, 

where  x  denotes   the  vector   (x„ ,x_ x  )'    and  V  is   the  dispersion 

—  12m 

matrix  given  by 


V  = 


fl       p    ...    p" 
p         1   ...   p 

p        p    ...    1 

Due  to  the  simplicity  of  the  dispersion  matrix,  the  determinant  |v| 
is  easily  evaluated  and  the  quadratic  form  x  V   x  has  a  simple 
scalar  representation.   In  fact,  it  will  be  shown  that 


(1.2.3)      |v|  =  (l-p)m"1[l+(m-l)p], 


(1.2.4) 


-1 


(l-p)[l+(m-l)P] 


"l+(m-2)p       -p 
-0       li-(m-2)p 

L   -p         -p 


-p 
-p 

...   l+(m-2)p_ 


Since  the  determinant  of  a  matrix  equals  the  product  of  its 
latent  roots,  (1.2.3)  can  be  proven  by  finding  the  m  latent  roots  of  V. 
If  X  is  a  latent  root  of  V  it  must  satisfy 

(V  -  \I)y  =  0  , 

where  0  represents  the  (mxl)  vector  of  zeros,  for  at  least  one  non-zero 
vector  y.   Letting  X  -    (1-p)  we  note  that  the  matrix  (V  -  \I)  contains 
only  one  distinct  element.   Therefore  the  rank  of  (V  -  \I)  is  one  and 
there  exist  (m-1)  non-zero  vectors  y  satisfying  (V  -  >.I)y  =  0.   Hence, 
(1-p)  is  a  (m-1) -fold  latent  root  of  V.   To  find  the  last  latent  root 
we  note  that  the  trace  of  a  matrix  equals  the  sum  of  its  latent  roots. 
Since  the  trace  of  V  equals  m, 


\     =   m  -  (m-1)  (1-p) 
m 

=  1  +  (m-l)p , 


and  (1.2.3)  is  proven. 


-1 


To  prove  (1.2.4)  denote  V   by  A  =  (a  _).   Then  A  must  satisfy 

op 


m 


£  "bj  'a  •  *• 


m 


E  v  .  a.a  =  0,   a^P 
J=l  aj   J0 


Substituting  for  v  .  in  the  above  equations,  we  have 
(1.2.5)     p   Z  a  R  +  aBg  =  1, 


p  =  ajP  +  (1^)aap  +  paPP^  °- 


Subtracting  these  last  two  equations,  we  have 


(1.2.6) 


a80  ~  aa9  -  1-p 


It  is  clear  from  the  preceding  equations  that  a  fl  is  independent  of 
Ot   and  0.   Hence,  equation  (1.2.5)  simplifies  to 

(1.2.7)      (m-l)pa^  +  a  g  =  1. 

Finally,  the  solutions  of  equations  (1.2.6)  and  (1.2.7)  are 


O0  "  d-p)[l+(m-l)p]  ' 


1+  (m-2)o 


Bp"  "  (l-p)[l+(m-l)p]  ' 


as  was  to  be  shown. 


-1 


Using  the  scalar  representations  of  |v|  and  x  V  x  in  the 
density  function,  (1.2.1)  may  be  written  as 


(1.2.8) 


,  m 

P    (p)  = 
r :  m      \r 


m      m-1         _  1 
(2n)  2  (1-p)   :  [l+(m-l)p] 


[r,m] 


exp 


m 


[l+(m-2)p]  E  x.  -  2p  2  E  x.x. 

i=l   X      i<J   1  J 

2(l-p)[l+(m-l)p] 


dx  .  .  .  dx 
1      m 


Unfortunately,  although  void  of  matrix  notation,  the  above  representation 

of  P    (p )  is  not  too  desirable  for  obtaining  numerical  results.  A   more 
r:m 

workable  form  of  P    (p)  is  needed. 

r :  m 


1.3  A  Transformation  Simplifying  P    (p ) 

r :  m 


Consider   the  m  variables   x',x'..,x'    defined  by 

12  m 

Xx  =   Y1   -    9Yc 

X2   =  Y2   "    9Yo 

X'    =  Y     -    6Y    , 
mm  o 

where   the    (m+1)    variables,    Y    ,Y„ ,...,Y    ,    have   a  multivariate  normal 

o   1      m 

distribution  with  mean  vector  zero  and  dispersion  matrix  I,  and  6  is 

an  arbitrary  constant.   It  follows  that  x'  x' x'  have  a  multivari- 

1   2      m 

ate  normal  distribution  with 

E(x')  =  E(Y.)  -  0E(Y  ) 
1         1  o 

=  0, 

Var  (x')  =  Var  (Y.)  +  9  Var  (Y  )  -  29  Cov  (Y. ,Y  ) 
l  x  o  l   o 

=  1  +  92, 

Cov  (x!,x'.)  =  Cov  (Y.-9Y  ,Y.-9Y  ) 
i   J  l   o'  j    o 

=  Cov  (Y.,Y.)  -  9  Cov  (Y. ,Y  ) 
i   J  i   o 

2 
-  9  Cov  (Y  ,Y.)  +  9   Cov  (Y  ,Y  ) 
o  j  o  o 

=  92. 
Therefore,  assuming  p  >  0,  if  we  define  9  by 

9  =  _e_ 


(i-p)* 

x'      x'         x' 

2     2-1                                   1          2              m 
9  (1+9  )     equals  p  and  the  variables  ,  ,  .  .  .  ,  

(l+92)?   (1+92)2       (l+92)^ 
have  the  same  distribution  as  X  ,X  , ...,X  defined  in  Section  (1.2). 

X     A  Til 

That  is, 


Var 


cov 


(i.e2)* 


1 


i+e 
=  l, 

x' 


2   Var    (Xi) 


u+eV      a+eV 


L_  cov    (x'   x') 

2  l      J 


1+9 


1+9 


=  P 


Hence,   we  have 


P        (D) 
r:m 


r)p<xi>0 Xr>0;   Xr+1<0 Xn,<0) 


\<i+eV 


x 


(mV 


—  >  o- 


X  X 

r+i    <o,..., — 2L_ 


2  £ 
(1+9  ) 


(1+92)* 


-© 


P(X'>0 x'>0;    x'       <0,...,x'<0) 

1  r  r+1  ra 


=    (!)  P(Y  -9Y    >0,...,Y   -9Y    >0; 
\r/  1        o  r        o 


Y      ,  -9Y    <C Y     -  9Y    <0) 

r+1  o  mo 


This   last    expression  can  be   simplified  by   noting   that   the  variables 

Y    ,Y„ , . . . ,Y      are  mutually    independent  with   common  density   and  distribu- 
o      1  m 

tion  functions 


i    2 

1        ~zY 

g(y)      =     e  J  ,      -CO<y<CD) 


-/2TI 


G(y)   =  f  g(t)dt, 


respectively.   Hence,  writing  P    (p )  conditional  on  Y  =  y  and  inte- 
f  j  r :  m  o 

grating  over  y,  we  have 

Vm^  »  (r)/"  P(Yl>eYo'---'Yr>9Yo;  Yr+l<9V 

_00 

..  .  ,Y  <6Y  |y  =y)  g(y)dy. 
m    o '  o 

Finally,  using  the  independence  and  identically  distributed  properties 

of  Y„ ,Y_ . . . . ,Y  ,  we  have 
12      m 

00 

(1.3.1)     P   (p)  =  (m)    f  [l-G(ey)]r[G(ey)]ra"r  g(y)dy. 

r:m      \r/  J 

—  CO 

Results  similar  to  (1.3.1)  have  been  given  by  Ruben  (1954),  Dunnett  and 
Sobel  (1955),  Moran  (1956),  and  Stuart  (1958). 

Although  the  expression  for  P   (p )  given  in  (1.3.1)  was  derived, 

r  :m 

assuming  p  >0,  Steck  and  Owen  (1962)  have  shown  that  it  also  holds  for 
p  <0  by  defining  G(6y)  in  the  complex  plane.   For  p  <0,  9  is  an 
imaginary  number  and  can  be  written  L®     where 


cp  =  — e- 


(1-p) 

Then  G(8y)  equals  G(-t.<?y)  and  is  defined  by  integrating  along  a  path 
in  the  complex  plane  parallel  to  the  x-axis  from  -«  +  i.cpy  to  + -£.cpy. 
That   is, 

2   2 

iCp   y        o  . 

(1.3.2)  Gt&py)   =   e  J       e"^Py  g(t)dt- 


The  proof  of  (1.3.1)  for  p  <0  consists  of  showing  that  the 
right-hand  sides  of  equations  (1.2.8)  and  (1.3.1)  are  identical. 
First  we  note  that 

,2„2   o 


1  -  GCi.9y)  =  1  -  e8     J 


•**»  g(t)dt 


i  2  2   °° 
l-/y  J   e-^  g(t)dt 


L.„2  2   » 


49  y   P   -ittpy 
+  e*     J   e     g(t)at 


,2  2" 


J9y   p   -i.tcpy   ...  .. 
!*•     J   e      g(t)dt. 


Then  using  (1.3.2)  and  writing  the  right-hand  side  of  equation  (1.3.1), 
say  R,  as  a  multiple  integral,  we  have 


00 

R  =  £)  J    [i-G(icpy)]r[Gacpy)]m~r  g(y)dy 


-©J. 


1-1 

J  <  r 


r   2   2  - 

-ff  y        -itpy     £     tj 

e  e  J=1         g(t   )...g(t   )dt    ...dt 

1  r        l  r 


m 


m-r  2   2      -Xpy        £        t . 


J— J  e 

e 

t.  <0 

j>r 

»     m2   2 

UJ     e 

-00 

J    • 

^[r,m] 

j-r+1       g(t     ,)...g(t   )dt     ....dt 


r+1' 


m        r+1  m 

g(y)dy 


m 


-l<py  £     t 
j  =  l 


g(t1)...g(t   )dt.. .dt 

1  m        1  m 


g(y)dy. 


10 


Since     p   > 


-1 
m-1    ' 

2        -mp 


r*D 


1-P 


m 


m-1+1 


=   1 


and   the  integral 


m 


-  y,*  ~*yA  *> 


g(y)dy 


converges.   Therefore,  interchanging  the  order  of  integration  is 
permitted  and  we  have 


m 


-C)J 

[r,m] 


-  -iy(cp  S  t  )  2  2 

J=1  J   i  "i(i-V)y 


P 

J  e 

_  —00 


-v/StT 


dy 


g(t)...g(t  )dt    .dt 
1       ml     m 


2  £ 
The  integral  in  brackets  multiplied  by  (1-ntp  )   is  the  characteristic 

function  of  a  normal  random  variable  with  zero  mean  and  variance 

2  -1 

(1-ntp  )   .   Hence, 


R  = 


m    2 
(9  E  t  ) 

m   2     1-1  J 

(r)  I   <*T)-*(lV>~*  e   J"1   J  e2(1-^>dtr 
[r,m] 


dt 


m 
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-1  2 

Substituting     -p(l-p)        for  9    , 


m        2 

m  -4     E     t 

J 


R  =    (™)     J      (2TT)"2(l^)"*[l+(n,-l)]"*   e        j=1 
[r,m] 


/  m  ,2 

p(Z     t 

e2[1+(m-1)p]      dt   ...dt    . 
1  m 


Finally,  making  the  transformation 


x . 
t.  =  — -2-r  ,   j  =  1,2,...  ,m, 

3     a^)s 


we  have 

m      m-1 


R=  Q     (2n)  2(1^3)     [l+(m-l)p]  * 


[r,m] 

m 


2 


[l+(m-2)p]  E  x.  -  2p  E  S  x  x  "\ 

-  ^ i^ }dx  .  .  .  dx 

2(l-p)[l+(m-l)p]        J  m 

and  the  proof  is  complete. 

The  results  of  this  section  for  positive  and  negative  rho  can 

be  summarized  in  the  following  lemma. 


Lemma  1 


03 

?r.m(p)  -  (")  /  [l-G(9y)'jr[G(ey)]m"r  g(y)dy 


where 


e 1 

(i-p)5 


12 


and  is  written  as  L®     when  p  <0  with  9  defined  as 


cp  =  ±!2j5 


(1-p)* 


The  functions     g(y)      and  G(9y)      are  defined  by 

2 


g(y)    =    —    e    "        ,        -oo<   y   <    co 


2TT 


Leo    *(t)dt,  P   >  0 


G(9y)      =  I 


^  L  *~^y  e(t)dt,        p  <  o 


^e 


1.4   Summary  of  the  Results  of  This  Dissertation 

In  the  next  chapter  Lemma  1  and  a  resuL.it  connecting  the  terms 

of  a  binomial  series  with  Tchebycheff-Hermite  amd   Legendre  polynomials 

are  used  to  obtain  a  finite  series  expression  fk*r  P   (p )  valid  for 

r  :m 

all  allowable  values  of  p.   For  p  <0  the  resudlts  reduce  to  a  work- 
able form  once  the  real  and  imaginary  parts  of   G(i.cpy)  are  isolated. 
In  Chapter  3  it  is  shown  how  the  results  of  Chapter  2  can  be 
programmed  to  obtain  numerical  results.   Since   P  (p )  is  simpler  than 

P    (p  )  for  computing  purposes,  a  result  expressing  P    (p }  as  a  sum 
r  :ia  r  :m 

involving  P .  (o  )  ,  .j  =  r,r+l,...,m,  is  proven.   "The  accuracy  of  the 
computed  results  is  verified  by  comparison  with  exact  results  for 
special  values  of  m  and  p  and  with  the  results  <*i>f  Ruben  and  Gupta. 
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An  application  of  the  results  of  the  first  three  chapters  is 
given  in  Chapter  4.   The  statistic  suggested  by  David  is  used  for  test- 
ing the  hypothesis  that  a  population  is  normally  distributed.   The  null 
distribution  is  discussed  and  numerical  resuJ ts  are  given  for  sample 
sizes  not  exceeding  22.   Approximations  to  the  null  distribution  are 
also  given. 

Finally  in  Chapter  5  we  discuss  three  alternative  methods  for 

computing  P  (p ) .   However,  none  of  these  methods  can  be  used  to  obtain 
m 

numerical  results  as  readily  as  the  method  discussed  in  Chapters  2  and  3. 


CHAPTER  2 

AN  EXPRESSION  FOR  P    (p )  INVOLVING  TCHEBYCHEFF-HERMITE 

r  :m 

AND  LEGENPRE  POLYNOMIALS 

2. 1  Definitions  and  Properties 

Let   c  (r,m)  denote  the  kth  order  Tchebychef f-Herir.ite  polynomial 
orthogonal  on  r  =  0,l,...,m.   Then   c  (r,m)  can  toe  written  (see,  for 
example,  Plackett,  Sec.  6.5),  as 


(2.1.1)      ck(p,m)  =T555T  E  (-1)   Vk_jy^  k  X  j  )    ' 


and  satisfies   the  following   three  properties: 


m 
(2.1.2)  E     c    (r,m)    =   0,  k=l,2,..., 

r=0     K 


m 
(2.1.3)  E     c.(r,m)c   (r,m)   =0,  j   /  k, 

r=0     J  k 


m  rvM2   /'m+k+1 

(2.i.«       I  oJ<r..> .:      u^- 


r=0 


/2k\ 


Also,  let  L  (t)  represent  the  kth  order  Legendre  polynomial 
in  t.   L  (t)  is  given  by  (see,  for  example,  Abraraowitz  and  Stegun, 

K 

k  2  ~h 

Chap.  22)  the  coefficient  of  s   in  the  expansion  of  (l-2ts+s  )  "    and 

can  be  computed  from  the  recurrence  relation 
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15 


L  (t)  =  1 

o 

(2.1.5)  L  (t)  =  t 

2k  — 1  k— 1 

Lk(t)  =^"t  •  Lk-l(t)  "IT  Lk-2(t>>      k*2- 

The  following  result  due  to  Saw  and  Chow  (1966)  connects  the 
terms  of  a  binomial  series  with  the  Tchebychef f-Hermite  and  Legendre 
polynomials.   For  any  p, 

m  2 

(2.1.6)  ^  U   P  (1"P)     Ck(r'm)  =  (m-k)!(2k):  \VP-»- 
r=0 

The  importance  of  this  result  to  the  next  section  justifies  the  inclu- 
sion of  the  following  proof. 

After  substituting   c  (r.m),  as  defined  by  (2.1.1),  into 

equation  (2.1.6)  and  multiplying  the  equation  by  — ,  the 

(kl) 

equation  to  be  verified  simplifies  to 

m   k 

V   V  f   11J  fm\  (  r  \   (Sk-JlUm-k+j):   r,,   ,m-r  ,T  ,_  ,, 

^  ,4  W  W  "    j:(k-j): p  (1"p)     =  ^V2p-1} 

Notice  that   I    J  equals  zero  unless   r  >  k  -  j.   Thus  letting  Q 
represent  the  left-hand  side  of  the  above  equation  and  changing  the 
order  of  summation,  we  have 

j=0  r=k-j 

Letting  r  =  r-k+j ,  the  sum  in  brackets  simplifies  to 


n-k+j 


Zm!  r  +k-j    ,„      .m-k+j-r 
p  J    (1-p) 

(m-k+j-r')  Ir'Uk-j)' 


r    =0 


m-k+j 


(k-j)pk_j     I       (m-^J)Pr/(l-p)ra-k+J-r' 


r    =0 


•  O  >-  • 


Therefore  Q  reduces   to 
k 

ji(k-j)!(k-j): 


Q .  l  (_uj  ^;,^>^,  Pk- 


Since  by  definition 


00 


k  r  2~i 

S     s  L    (2p-l)    =      l-2s(2p-l)+s 
k=0  *  L  J 
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it  remains   to   show  that 

1      °°       k         r  2~!~£ 

4-     E     S  Q  =      l-2s(2p-l)+s 

m-    k=0  L  J 

Making   the   change  of  variable      £  =   k   -  j      the  left-hand   side  of   the 
above   equation   becomes 


4=0  j=0 


,vj    (2i+j)I    nO  +  i 

j :  T.  V.  p 


2  -4 
However  it  can  be  shown  that  [l-2s(2p-l)+s  ]  s      also  reduces  to  the 

above  sum.   We  have 
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9       _1 

[l-2s(2p-l)+s"]      s   =    [i-4sp+2s+s    ]      2 


2.  ~h 


=    [(1+s)2   -  4sp]  ~7 


=    (1+s)  1 


'E1-  tttJ 


-* 


(1+s) 


=    (1+s)  E     (-1)      I     /=  J 57 

4=0  W  '    (l+s)2X 


/  _i\                                                            jg      (2i)  ! 
But      (     «     )  may  be  written  as      (-1)      ^— j 


(4!)    4 


since 


f-i\  _     V      2A      2/---V  2 


I 


24-1 
2 


x  y 


i: 


=  (-D 


4   1-3    .  .       (  21-1) 

2£  4: 


X  (24)! 

=    (-1)      -j —^ 


=    (-D 


2      X:    2.4    ...    24 
4      (21).' 

4IaT7  ' 


2,  -4 

Substituting   the  above  result   into  the  expression  for    [l-2s(2p-l)+s   J         , 

we  have 


[l-2s(2p-l)+s    ] 


21  -J  _      (2i)I     .      .4  1 


(sp) 


4=0      (X!)"  (ln-S) 


2TTi 


S     1E4):     (sp)£     E(-1)J    /24k,  N     J 

juo  (4.*)  j=o  v   J    y 


v       v     /■    i^J    <24+j)!      4  cj+4 

,En  .En  (_1)    i:4!4:    p    s      - 

4=0  j=0 


which   completes   the  proof. 
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2.2   The  Fundamental  Result 


In  this  section  we  use  the  results  of  the  last  section  and 

Lemma  1,  page  11,  to  obtain  a  finite  series  expression  for  P    (p ) . 
'   °  r :  m 

First  we  let  p  =  1  -  G(8y)   in  equation  (2.1.6).   Then  we  have 


m   ,  . 

y     !m\r.    „,  «__*  -irr  „,  „_.x.,m-r 


r=0  ^r 


j[l-G(9y)]  [GOy)]"1   c^r.m) 


m:(k:>2-rLv[l-2G(ey)] 


(m-k)!(2k).'   k' 


Multiplication  of  the  above  equation  by   g(y)   and  integration  with 
respect  to  y  gives 


00   m   ,  . 
(2.2.1)     J*    E  ^J[l-G(9y)]r[G(9y)]m"r  ck<rfm)  g(y) 


m 

E 
-  »  r=0 


2     °° 
-   f   L.  [l-2G(6y)]  g(y)  dy 


(m-k) ! (2k) 


_  CO 


Now  taking  the  left-hand  integral  inside  the  summation  and  defining 
Jk(p)  by 


(2.2.2)     Jk(p)  =J    Lk[l-2G(Gy)]g(y)  dy; 


2    2  -1 
where  p  -    6  (l+8~)   ,  equation  (2.2.1)  becomes 


m  \ 

E  ck(r,m)  /   ^)[l-G(6y)]r[G(6y)]m  rg(y)  dy 


m 

E 

k=0 

2 


m!(k!)" 

r  J,  (p) 


(m-k)! (2k)!   k 


19 


Finally,  applying  Lemma  1  to  the  left-hand  side  of  the  above  equation, 
we  have 


(2.2.3)     ^  ck(r,m)  P^o)  =  j^^yr   ^(p) 


An  alternate  expression  for  P  _  (p )  can  be  obtained  by  noting 

that  for  fixed  m  the  set  of  points  P   (p ) ,  r  =  0,1,..., m,  lie  on 

r:m 

a  polynomial  of  degree  at  most  m.   Hence,  for  some  constants,  say 

e  , e„ ..... e  ,  we  can  write 
o   1      m 

m 

(2.2.4)  P   (o)  =  E  e  c  (r,m). 

Multiplying  (2.2.4)  by  c  (r,m)  and  summing  over  r  yiel 

K 

m  mm 

(2.2.5)  E  c.(r,m)  P   (p )  =  Z   E  e  .c .  (r  ,m)  c,  (r  ,m)  . 
r=0  k       r:m      r=0  j=0  J  J      k 

Using  the  properties  (2.1.3)  and  (2.1.4)  of  c  (r,m),  the  right-hand 
side  above  reduces  to 


ds 


m   2 

6k   E   ck(r,m) 
r=0 


\ 72ky 


Hence,  the  constants  e  , e   ...,e   can  be  determined  by  equating  the 

o      l  m  ° 

right-hand   sides   of   equations    (2.2.3)  and  (2.2.5).       That   is, 

2   f'm+k+l\ 
(k,)      Vjk+jJ  m.'(k:)2 

k  72k\  -   (m-k) !  (2k)  :   Jk(p)- 


(f) 


20 


After  slight  simplification,  we  have 


f  m!  (2k+l):  1  T  ,  . 

=  L — * JJk(P> 

'-(kl)  (m+k+1)'.-1 


ek  =  I o !J,.(P).         0,1,..., m 


Letting  h  (m)  denote  the  constant  in  brackets  and  substituting  into 
equation  (2.2.4),  we  have 


m 

(2.2.6)     P   (p)  =  Z     c.(r,m)  b,  (m)  J,  (p ) . 

r:m      ,  .   k       k     k 
k=0 

Before  investigating  the  integral   J  (p ) ,  we  should  comment  on 

the  utility  of  the  expression  for  P    (p )  given  by  (2.2.6).   Most 

r  :m 

important,  by  defining  J  (p )  appropriately,  the  expression  is  valid 
for  both  positive  and  negative  values  of  rho.   Next,  the  integral  J  (p ) 
does  not  depend  on  r  or  m.   Hence,  for  a  given  value  of  p,  only  one 
set  of  values,   J  (p ) ,  k  =  0,1,...,  is  needed.   Also,  as  will  be  shown 
in  the  next  section,  J  (p )  =  0  for  odd  k,  thus  decreasing  the  number 
of  terms  in  the  series  by  one-half.   Furthermore,  for  large  m,  the 
series  may  be  truncated  without  serious  effect,  since  the  factor 
c  (r,m)  b  (m)  approaches  zero  as  k  increases. 


2.3  The  Integral  J  (p ) 

K 


Consider  the  integral 


Jk(p)  =  /   Lk[l-2G(9y)]  g(y)  dy 


defined  for  k  =  0,1,...   in  the  last  section.   Using  the  recurrence 

relation  (2.1.5),  it  can  be  seen  that  the  Legendre  polynomial  L  (t)  is 

k 
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an  even  or  odd  function  in  t,  depending  on  whether  k  is  even  or  odd, 
respectively.   Also,  when  rho  is  positive,  G(8y)  is  the  normal  distri- 
bution function  which  implies   l-2G(9y)   is  an  odd  function  in  y. 
Therefore,  L  [l-2G(8y)]  is  an  odd  function  in  y  when  k  is  odd  and  an 
even  function  in  y  when  k  is  even.   Since  g(y)  ,  the  normal  density 
function,  is  an  even  function  in  y,  it  follows  that  J  (p )  ,  rho  posi- 
tive, equals  zero  when  k  is  odd,  and  for  even  k 


Jk(p)  =  2  J  L  [2G(6y)-l]  g(y)  dy. 
o 


After  making  the  transformation  u  =  6y,  J,  (p )  becomes 

K 

2 

oo  ,  u 

2  r  i    ~ft  oz 

Jk(p)  =  q  j  Lk[2G(u)-l]  -=-   e    -      du  . 

V2n 
o 

Finally,  letting   h  =  —g   and  writing  the  integration  with  respect  to 
the  distribution  function  G(u) ,  we  have 

1  1  2  h  , 

J  (p)  =  2h*   L  [2G(u)-l]  I  e  dG(u)  , 


e2        i 

where  p  =  ^  =  TTT  and  0  <  p  <  1  imply  that  h  >  0. 


Next  consider  the  integral   Jk(P)  when  rho  is  negative.   Now 

i  i 

8  =  -r  is  imaginary  and  is  written  as   8  =  &p   where  cd  =   ~^ — r  . 

<l-p>*  (1-p)* 

Hence,  the  function   G(9y)   appearing  in  J  (p )  is  complex  and  in  order 
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to  simplify  J  (p )  the  real  and  imaginary  parts  of  GCLyy)    must  be 
isolated.   Denote  these  real  and  imaginary  parts  by  Qf(cpy)  and  |3(cpy), 
respectively.   Then  GCi9y)  and  its  conjugate  can  be  written 

Gulcpy)  =  a(?y)  +  i.P(cpy)  , 
(2.3.1) 

G(-^9y)  =  a(9y)  -ISiyy). 

Using  the  definition  of  GCi.9y)  given  in  equation  (1.3.2),  we  have 
2a(cpy)  =  e*9^  {J  e-^y  g(t)dt  +  J  /*"*  g(t)dt} 

—  00  —00 


^2y2  J   e-^tcpy  g(t)  ^ 


2  2 

e^ 

_  00 


,22      i  2  2 

=  e?9  y  .  e-*P  y 


=  l. 

Therefore, 

or(cpy)  =  ^  • 

Subtracting  equations  (2.3.1)  gives 


2  2    ° 


2ig(cpy)  =  e^y  {J  .*^g(t)dt-J  eXCPyt  g(t)dt} 


i  2  2 

=  e^9  y 


'  (e^  -  e^)  g(t)dt. 


Since  e     -  e     =  2j_  sin  (tcpy)  ,  we  have 


i  2  2  * 
(2.3.2)     0(cpy)  =  e*9  y  J*  sin  (tcpy)  g(t)dt 
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3  (cpy)  can  be  further  simplified  through  integration  by  parts 
and  differentiation  of  (2.3.2).   First,  integrating  by  parts,  we  have 

j.  2  2  -  -iro 

(2.3.3)     3  (cpy)  =  e2^  y   |-  —  cos  (tcpy)g(t)J 

00  . 

-  £-  J   t  cos  (tcpy)g(t)dt| 


?y  o 


±  .    _1_  e*pV  _  ±  &V  ft   cos  (ttpy)g(t)dt. 


Next,  differentiating   S(cpy)  with  respect  to  cpy,  yields 


i  2  2   °= 

(2.3.4)     *,*?  =  csy  e^  y   f  sin  (ttpy)g(t)dt 
d(Vy)  o 


2  2° 


+  e*^  y  J  t  cos  (tcpy)g(t)dt, 


where  differentiation  was  permitted  inside  the  integral,  since 

jt  cos  (tcpy)g(t)|  <  tg(t)   which  is  integrable.   Combining  equations 

(2.3.2),  (2.3.3),  and  (2.3.4),  we  have 

«<&>  =  cpy  0(cpy)  4  -i_  e^V  -  97  e(W) 
d(Cpy)  v^ 

i  2  2 

=  _L.  e^P  y 

V^TT 
It  follows,  since   0(0)  =  0,  that 


<?y    x      it2 

(2.3.4)     8  (cpy)  m   j    e2    dt. 

o   72tt 


Hence  the  complex  function  G(-iipy)  can  be  written  as 
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cpy  i    2 

GOJPy)    =  -  +  jl  J         eE        dt 

o       v^2tt 


=  |  +  i.0(cpy) 


Returning   to   the   integral      J^(p)    we   can  now  write 


Jfc(p)    =  I     Lk[-2XS(cpy)]    g(y)    dy. 

_00 

Although  the  Legendre  polynomial  has  an  imaginary  component,  its 
definition  and  recurrence  relation  still  hold.   In  fact, 
L ,  [ -2ij3  (cpy )  ]   is  an  even  function  in  y  when  k  is  even  and  an  odd 
function  in  y  when  k  is  odd.   Therefore,  as  in  the  case  when  rho  is 
positive,  J  (p )  equals  zero  for  odd  k  and  for  k  even 


J  (p)  =   2  J     L  [2Z0(9y)]  g(y)  dy. 


k 
o 


After  making  the  transformations  u  =  cpy  and  h  =  —  =  -  — 

8      <p" 


J,  (p )  becomes 
k 


-h-1 


i  l  r  jl.2-1 

Jk(p)  =  2(-h)*  J*   hk\2l   P(u)]  I  e-6"        dG(u), 
8 


where  p  =  ; — -  and   -  <  p  <  0   imply  that  h<-m,  m  ^  2. 

h+1         m-1 

The  following  lemma  summarizes  the  results  of  this  chapter. 


Lemma  2 
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P    (p) 
r:m 


m 

I    c  (r,m)  b  (to)  Jk<P>, 
k=0 

k  even 


where 


V-B)  ■  w  j0  (-1)J 


U  (T)  C 


m-k+ j  \ 


b.  (m) 

k 


m!  (2k+l)I 
(kl )2(m+k+l): 


and  for  even  k,  J.  (p  )  is  the  integral 
k 


with 


i   1         j-   i  2-,|h|-l 
Jfc(p)  =  2Jh|2  /  Lfc[h(u)]  \j~  *   X  dG(u) 


P  = 


h+1 


and 


h(u)  =  / 


2G(u)  -  1 


%  s  .^ 


h  >  0 


v.V21t   o 


ed        dt     h  <  -m,  m  ^  2 


L  (t)  is  defined  by  the  recurrence  relation 

L  (t)  =  1 
o 


Lx(t)  =  t 


2k  — 1  k— 1 

Lk(t)  =  -ir  t  \-i  (t)  -nrLk-2(t>-   ks2- 


CHAPTER  3 


NUMERICAL  RESULTS 


3. 1  Exact  Results 


In  general  the  value  of  the  integral  J  (p)  can  only  be 

K 

approximated  so  that  exact  results  for  P    (p )  are  not  available. 

r:m 

However,  exact  values  of  J  (p )  and  J„(p)  can  be  found.   Then,  since 

o         Z 

J  (p  )  is  independent  of  m  and  r,  P    (p),  r  =  0,1,2,  and  P    (p )  , 

K  r  *.  <i  r ;  o 

r  =  0,1,2,3,  can  be  determined. 

Since  L  (t)  =  1,  we  have  from  Lemma  2,  page  25, 
o 

i    i  i-  i  2-j  |  h  |  -l 

JQ(p)  =  2|h|?  J*  Lo[h(u)]|e"cU         dG(u) 

|h|   2 

r       h  [■    2  u 

=  J    - — <—   e         du 

=  1. 

J  (p )  can  be  determined  indirectly  by  first  finding  P  (p ) .   Letting 
m  =  r  =  2  in  Lemma  2 , 

V        (p)  =  P  (p)  =  c  (2,2)  b  (2)  J  (p) 
&'.£.  £  o  o  o 


+   c2(2,2)  b2(2)  J2(p), 


so  that 


P.(p)  -  c  (2,2)  b  (2) 
(3.1.1)     J0(p) 


2V'  c_(2,2)  b„(2) 
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The  value  of  P„(p)  can  be  found  in  closed  form  by  integrating  the 

original  expression  for  P  (p )  given  in  equation  (1.2.8).   With 

m 


m  =  r  =  2 ,  we  have 


P2(p)  - 


(2^)  *  (1-pV2  e 


.  2(1^2)    -v* 


o  o 


Making  the  transformation 
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xrpx2 


d-p  )■ 


U2  =  X2  ' 


it  follows  that 


2  & 
x1   =  (1-p  )*  ux  +  pu2 


X2  =  V 


and  the  Jacobian  of  the  transformation,  J(x  ,x  —  u  ,u  ) ,  is 


J(X1'X2  "*  U1'U2}  = 


2  4 
(1-P  )     P 


-  (I-.2)* 


Therefore, 


P2(P) 


2*  6  dU2?dUl 
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Finally,    making   the  polar   transformation, 


u     -   r   cos    0 


u     =   r   sin   6, 


we  have 


°°    r      cos      (-p)  .  1 


=  —  cos      (-p) 

1    /tt  .    -1     , 

=  2Tl2+Sin       P) 

1        1  -1 

=    4  +    2^  Sin       P  • 

Next  we  need  the  values  of   c,  (r,m)  and   b.  (m)  for  k  =  0   and 

k  k 

k  s  2.   For  k  =  0,  c  (r,m)  =  1   and  for  k  =  2, 

o 

c  (r,m)  =  r(r-l)  -  r(m-l)  +  -  m(m-l) . 

12      1 
Thus  for  m  =  2,  c  (r,m)  equais  -^,    -   — ,  and  —  for   r  =  0,1,  and  2, 

respectively,  and  for  m  =  3,  c  (r,m)  equals  1,  -1,  -1,  and  1  for 

r  =  0,1,2,  and  3,  respectively.   The  constants  b  (m)  and  b  (m)  are 

given  by 

o      m+1 

b0(m; 


2V      (m+1) (m+2) (m+3) 


Substituting  m  =  2   and  m  =  3  gives 


bo(2)  =  \  ■         bo(3)  -  i  • 

b2(2)  =  |  ,        b  (3)  =  |  . 
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Now  using  (3.1.1)  we  find  that 

J2(p)  =  -  sin   P  ~  2  ' 
Finally,  substituting  the  above  results  into  Lemma  2,  we  have,  for  m=2, 

P0:2(p)  =1+   2^Sin_lp 

1   1   .  -1 

Pl:2(p)  =  2  "^Sin   P 

P2:2(P)  =\+   2^S1I1"lp* 


and  for  m  =  3 , 


1    3.-1 

P    (o )  =  —  +  —  sin   p 
0:3VM'    8   4tt       m 


P1:3(P)  =1"  J?Sin~lp 


P2:3(P)  =f-  4^Sin_lp 


,      ^      1  3  _1 

P3:3(P)  =  8+  4^Sin   p- 


Notice  in  the  above  results  that 


P„   (p)  =  Pm  v.  (P>> 
r: m       m-r :m 


and 

m 

Z  P    (p)  =  1. 
r:  m 
r=0 

These  properties  also  hold  in  general.   The  first  result  follows  imme- 
diately by  letting  r  =  m-r  in  Lemma  1,  page  11,  and  noting  that 

G(9y)  =  l-G(-9y)  for  both  positive  and  negative  rho.   Lemma  2  is  used 

m 
in  showing  the  second  property.   Since  S   c  (r,m)=0  for  k=l,2, 

-H f\  K 

we  have 


r=0  k 
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m  mm 

E  P   (p)  =  E   E  b,  (m)  c  (r,m)  J  (p ) 


r : m        _  .  ~ 
r=0  r=0  k=0 


m 


=  E  b  (m)  c  (r,m)  J  (p ) 
,.oo       o 
r=0 


m   1 

=  E 


r:m 
From  Lemma  1, 


.  m+1 
r=0 


P  _(p)  can  also  be  computed  exactly  in  the  case  when  p  =  — 


CO 

P   4)  =  (m)  I      [1-G(y)]r  [G(y)]m_r  g(y)dy, 
r:m  2     \r/  "J 

—  00 

since   9-1  when  P  =  o  "   ^^ 

Q  [1-G(y)]r  [G(y)]m"rg(y) 

where  the  term  inside  the  braces  is  the  density  on  the  (m-r+l)st 
order  statistic  from  a  normal  random  sample  of  size  (m+1). 
Therefore, 

P   4)  =  rrr  .     r  =  0,1 m. 

r:m  2    m+1 

This  last  result  implies  that   J  (^-)  =  0,  k^2.   From  Lemma  2, 

1     m  1 

P   (±)  =  E  b  (m)  c,  (r,m)  J  C^) 

r:m  2    ,  _  k     k       k  2 
k=0 

k  even 

m 

=  — 7  +   E    b,  (m)  c  (r,m)  J  (-) . 
m+1     „    k     k       k  2 

k  even 
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Since  the  last  sum  equals  zero  for  all  r  and  ra  and  since  b  (m)    c  (r,m) 
does  not  equal  zero,  J.Ct,)  =  °- 


3.2  Evaluating  the  Integral  J,(p  ) 


The  Case  When  Rho  Is  Positive 

Recall  from  Lemma  2,  with  p  -    (h+1)   >   0,  that 

f*    A-  r    A   2— ;h-l 

Jfe(p)  =  2  J   h*  Lk  [2G(u)-l]    e  "  U        dG(u)  . 

Before  attempting  numerical  integration,  the  integrand  should  first  be 
investigated  for  different  values  of  k  and  p  .   For  given  values  of 
G(u) ,  the  Legendre  polynomials  can  be  evaluated  from  the  recurrence 
relation  (2.1.5).   Also,  the  value  of  u  corresponding  to  G(u)  is  given 

to  eight  decimal  places  in  The  Kelley  Statistical  Tables.   For  each  of 

11         9 
the  cases  p  =  —  ,  —  ,  and  —  ,  the  integrand  is  plotted  in 

Figures  1,  2,  and  3,  respectively,  for  both  k  =  2   and  k=10. 

Since  numerical  integration  is  most  accurate  when  the  function 
being  integrated  is  well  behaved,  we  should  expect  good  results  for 
p<  —  ,  with  the  accuracy  increasing  as  rho  decreases.   Also,  since  k 
is  the  order  of  the  polynomial  being  integrated,  the  better  results 
should  occur  when  k  is  small.   When  p  >—  ,  the  integrand  approaches 
infinity  as   G(u)  approaches  one  as  can  be  seen  from  Figure  3.   In  this 
case,  it  is  not  likely  that  numerical  integration  will  give  accurate 
results. 

Using  Simpson's  rule  with  200  intervals,  the  interval  width 
equals  —-  -    .0025,  with  G(u)  taking  the  values  .5  +  (.0025)j, 
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1.0 1-  3Lk[2G(u)  -1]  h 


r  _4 

e  a 


2-i  8 


J 


1.0 


G(u) 


i         r  ..iu  1 

Figure  1.   The  function  hs   I.  [2G(u)  - 1]   e  2 
p  =  1/10. 


2-ih-l 


33 


l.Oj.  Lfe[2G(u)  -1] 


5+- 


-4—  G(u) 


-1.0- 


1:  =  2 


k  ^  10 


-1.5 


l  2-,  h-1 


Figure  2.   The  function  hs  L,  [2G(u)  - 1] 


V 


[•■*ri 


p  =  1/2. 
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2-1-8/9 


1.0- 


-1.0  + 


i  r-_  i  2-,h-l 

Figure  3.   The  function  h2  I.  [2G(u)  -  1]   e  ^U 

p  =  9/10. 
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j  = 0 ,1, . . .  ,200.   Then  J.  (p )  can  be  approximated  by 

i  0025  20°  r-£u--]h-l 

J.  (p)  -  2h*  ^^  S   c.  L.(.005j)  |  e    J 

where 

c  =  c    =1 
o    200 


J 


4   ,  j  odd, 
c  . 

2      ,   j  even, 


and  u.  is  the  value  of  u  satisfying   G(u)  =  .5  +  (.0025)j. 

The  Case  When  Rho  Is  Negative 

With  p  <0,  the  Legendre  polynomial  has  the  imaginary  argument 
■if  (u)  ,  where  f  (u)  is  given  by 

f(u)  =  J-     f   e*1   dt. 

<n/2TT   O 

For  given  values   of  u,    the  function  f (u)    can  be  evaluated  quite  rapidly 
by   first   expanding   the   integrand   in  a  Maclaurin  series   and  integrating 
term  by   term.      That   is, 

v^n     o     j=0  2Jj 


—    J       £     ^-  dt 

^        -;_n    o-J  -j  ' 


2  ™      ^U   t2j 

-       E    j      ~—  dt 

■v/2TT     j=0   0      2 


j,. 


-    2       «*w 


J=0V5n     2J(2j  +  l)j: 
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Letting  f  .  (u)  denote  the  (j-«-l)st  term  of  the  sum,  we  have  the  recurrence 
relation 

f  (u)  =  —  u 

2 

f  .  (u)  =  ~  •  .  roJ"  \  :  f  .  .  (u)  ,      3  =  1,2,..., 
j       2    J(2j  +  1)   j-1 

which  can  be  easily  programmed. 

Although  the  Legendre  polynomial  has  the  imaginary  argument, 
Zf (u) ,  for  even  values  of  k  it  is  a  real  function.   Hence,  for  comput- 
ing purposes,  we  can  avoid  complex  numbers  by  defining  the  function 

*  . 
L  (-it)  by  the  recurrence  relation, 

K 

L*at)  =  i 

o 

(3.2.1)     L*(it)  =  -t 

*  .  k  2k-1        *  k— 1   * 

LkUt)  =  <"l)   -k-  *  •  Lk-Iat)  "  —  \-2Ut)>         k^2- 

Then  the  function  L,  (it)  can  be  determined  by 

k 

*  . 
L  (it)       ,   k  even 

Lkat)  =    *  . 

-iL,  (it)        k  odd. 
k 

The  above  relation  can  be  verified  by  substituting  into  the  recurrence 
relation  (3.2.1)  and  comparing  the  results  with  the  recurrence  relation 
for  L  (t)  given  by  (2.1.5). 

Once  the  Legendre  function,  L  [if(u)],  is  evaluated,  we  can 

i       ^.;u2vh-i 

examine  the  integrand,  (-h)8  L  [if(u)]  [  e  K  ,  of  J  (p ) .   The 

integrand  is  plotted  in  Figure  4   for  p  =  -—  ,  k  =  2,  and  for 
p  =  -  —  ,  k  =  2  and  k=10.   Notice  that  the  scale  of  the  graph  differs 
from  those  in  Figures  1,  2,  and  3,  since, for  large  k, |j  (p)|  is  quite 
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5  [.  (-h) s  LkUf(u)]  I 


r.-*u 


2-,-h-l 


-5  - 


-10  ■. 


-15 


-20 


J 


p  =  -  1/2,  k  =  2 

p  =  -  1/10,   k  =  2 

p  =  -  1/10,   k  =  10 


Figure  4.   The  function  (■ 

p  =  -  1/2,  p  = 


-h)s  I^DifCu)]  [_e"su 


2-,-h-l 


1/10. 
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large.   Also,  unlike  the  case  when  rho  is  positive,  the  functions  do 
not  cross  the  G(u)  axis.   These  differences,  together  with  the  smooth- 
ness of  the  curves,  should  make  numerical  integration  even  more  accu- 
rate for  p  <  0. 


3.3   Computing  P  (p ) 
m 


Once  the  J  (p),  k  =  0,2 m,  are  computed,  P    (p )  is  deter- 

k  r:m 

mined,  for  r  =  0,1,..., m,  by  the  expression  giveni  in  Lemma  2.   That  is, 

m 

(3.3.1)  P   (p)  =  Z  b,  (m)  c  (r,m)  J,  (p )  . 

r:m      ,  _  k     k       k 
k=0 

k  even 

Unfortunately,  since  the  Tchebychef f-Hermite  polynomial,   c  (r.m), 

cannot  be  expressed  in  a  recurrence  relation,  it  is  difficult  to  compute 

P    (p )  using  the  above  expression,  except  for  simall  values  of  m. 
r:  m 

However,  with   r  =  m,  c  (r,m)  simplifies  to,  say  c  (m)  , 

K  K 

(3.3.2)  c,  (m)  =   m:(k!) 


k      (m-k): (2k)!  ' 

Furthermore,  the  probability   P    (p )  can  be  written  as  a  linear  combi- 

r:  m 

nation  of  P  (p  )  ,  k  =  r,r+l,  .  .  .  ,m.   That  is , 

<3-3->     w>  -  C)  "v-»J  (7)  w»  • 

(ITie  proofs  of  (3.3.2)  and  (3.3.3)  are  given  at  tlie  end  of  this  section.) 

Hence,  we  use  (3.3.1)  only  in  the  special  case  wh.en  r  -  m.   Relation 

(3.3.3)  then  can  be  used  for  computing  P    (p )  wihen  r  ^  m. 

r :  m 
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A  further  simplification  in  the  computation  of  P  (p )  can  be 
made  by  combining  the  constants  b  (m)  and  c  (m) .   Letting  d  (m)  denote 
the  product ,  we  have 

d  (.m)  =  b  (m)  c  (m) 

m:C2k+l):       mldcl)2 


Ck!)  (m+k+1):    (m-k)  '.  C2k)  '. 


Cm:)2C2k+l) 


(m+k+1): (m-k) 


The  constant  d  (m)  can  be  computed  for  even  values  of  k  by  the  recur- 
rence  relation, 


(3.3.4) 


d   (m)    =  J_ 
o  m+1 


2k+l        m-k+2        m-k+1    .        ,    _  .         „    . 

di  (m)  =  tt, — o  •  — ; — r  •  ; —  d,  r,(nO  ,        k  =  2,4 m. 

k      2k-3   m+k+1    m+k   k-2    ' 


Finally,  combining  the  above  results,  we  have  the  following  computing 
formula  for  P  (p  )  : 


m 


171 


(3.3.5)     P  (p)  =   E  d,  (m)  J.(p  ) . 
m      ,  _   k     k 
k=0 

k  even 


The  proofs  of  (3.3.2)  and  (3.3.3)  follow.   To  prove  (3.3.2), 

we  first  let   r  =  m   in  the  definition  of   c,  (r,m\)  given  by  (2.2.1). 

k 

Then,  writing  combinations  as  factorials  and  cancelling  like  terms, 
we  have 


(kI)    .      m:      .    e    (-dj         (2k-j): 
(2k):      (m-k):      J=o  tik-j)!(k-j)!j! 
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Thus,  we  must  show  that  the  summation  above  equals  one.   This  summation 
can  be  written  as 

I  ^  (V)  (£)  ■ 

Then,  letting  I  -   k  -  j  ,  the  result  to  be  proved  becomes 


k 


L)*"       )  (*)  =  1. 


4=0 
Next,  we  introduce  the  negative  binomial  and  binomial  identities  given 

by 

(l+a)      r=0 


«.,.,,    (i+d*-Xg)© 


Multiplying  the  left-hand  sides  of  (3.3.6)  and  (3. 3. 7),  gives 


(1+a)k        1  1  I     (-I)2  a1 


k     k+1    k         k 
a  (l+a)      a  (l+a)    a   fcO 

Equating  this  product  to  the  product  of  the  right-hand  sides  of  (3.3.6) 
and  (3.3.7),  we  have 

-r     S   (-1)   a  =  E   S  (-1)        L  a 
ak  4=0  r=0  4=0      V  r  '    W 

o  k 

Finally,  equating  coefficients  of  a   and  dividing  by  (-1)  ,  we  have 


k 

E 


I   «-»**  (?)©-!. 


as  was  to  be  shown. 
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In  proving    (3.3.3),    we  begin  with    the  basic   definition  of 


P        (p).      That    is, 
r :  m 


P         (p)    =    (m)  p(x„  >0,...  ,X    >0;    X         <0 X    <o) 

r:m  \r/\l  r  r+1  m        / 


Next,    let  A.    be   the  event    [X.>0]    and   denote   its   complement   by  A.. 

Also,    define   I      as   the  intersection     AA-...A,  .      Then     P        (p )    can 
k  12k  r:m 


be  written    as 


p      (p)  =  (m)  p(a,...a  a    ....a) 

r:m  \r/      \  1  r   r+1  m/ 


=    (m)p(l     -A   ...A  A       ....I  ) 

\r/\r        1  rr+1  m/ 

=  (:)[p(v-Kv--AA+i--\)J- 

Applying  de  Morgan's  rules,  and  since  P(I  )  =  P  (p ) ,  we  have 

p   (p)  =  (m)  j  P  (p)  -pU...a      U  A  .)\ 
r:m  r     \r/  L  r      \  1    r  .    r+j/J 

\J  *~  " 

=  (;)[pr(p)-pCu/1...ArAr+.)] . 

Finally,  using  the  formula  for  the  probability  of  a  union,  we  have  the 
desired  result.   That  is. 


p) 

j-m-r 

-LE 

p(a1.  . 

.  .A  A          ) 

r  r+j/ 

E    E    p(a 

i    <  j 
Jl     J2 

1...A] 

A          A 
"  r+31  r+J2 

+  E       E       E     p(a  . 
31     J2     J3 


.A  A       .A       .A 
r  r+Ji  r+j2  r+jg, 
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f(V...a  a  ,...aJ)  j} 

\  1    r  r+1    m/JJ 


-...«-!)"- 


■C)[*,«-W'wW«w 


3. 4  Accuracy  of  the  Results 


The  integrals  J  (p )  ,  the  constants  d  (m)  ,  and  the  probabilities 

K  K 

P  (p )  were  evaluated  with  double-precision  accuracy   using  an  IBM  model 
m 

360  computer.   The  computed  values  of  J  (p )  were  checked  against  the 
exact  values  for   k  =  0  and  k  =  2  and  for  p  =  —  .   For  k  =  0  and  k  =  2  , 
the  results  were  accurate  to  at  least  seven  significant  digits  for 
[p I  <  — ,  and  to  five  significant  digits  for   |p  |  =  —  .   J  (— )  was 
computed  accurately  to  six  significant  digits  for   k<  14  and  to  five 
significant  digits  for   16<  k  <  22.   Hence,  for  Jp | <  —  ,  we  would 
expect   J  (p )  ,  k>2,  to  be  accurate  to  at  least  the  fifth  significant 
digit.   The  computed  values  of  J  (p ) ,  k<  22,  are  given  in  Appendix  1 
for  p  =  -  ,  p  =  2(1)25,  and  for  p  =  -  -,    p  =  3(1)26. 

As  expected,  with  p  >  — ,  accurate  values  of  J  (p )  were  not 
obtained  using  the  method  of  quadrature  described  in  Section  3.2. 
Further  investigation  has  to  be  made  in  order  to  find  a  means  of 
evaluating  J  (p )  accurately  for  p  >  —  . 

K  £t 
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The  constants  d ,(m)  were  evaluated  exactly  using  the  recur- 
rence relation  (3.3.4).   They  are  tabulated  in  Appendix  2  for 
m<  25,  k  s  0,2 m. 

Finally,  the  probabilities   P  (p )  were  evaluated  using  the 

m 

formula  (3.3.5).   Results  for  p  >0  were  compared  with  Ruben's  tables 
and  were  found  accurate  to  at  least  five  decimal  places.   For  p  <0, 
Steck's  relation, 

Vp)  =  2  "  *     W  Pk  \l+(k-2)pj  Pm-k  VTTkp/  ' 
k  even 


was  used  in  making  comparisons.    Again,  the  computed  values  of  P  (p ) 

were  accurate  to  at  least  five  decimal  places.   P  (p )  is  tabulated  in 

m 

Appendix  3  for  p  =  -,   p  =  2(1)  (25)  and  2  <  m<  22,  and  for 

p  _ ,  p  =  2(1)21,  and  2  <  m  <  p. 

p  —   — 


CHAPTER  4 


APPLICATION:   A  TEST  FOR  NORMALITY 


4. 1   Introduction 


Consider  using  a  random  sample,  Y„,Y_,...,Y  ,  from  a  continuous 

12      n 

distribution  F  to  test  the  hypotheses: 

H  :   F  is  a  normal  distribution 
o 

H  :   F  is  a  skewed  distribution, 
a 

As  David  (1962)  suggested,  one  might  consider  the  number  of  observa- 
tions larger  than  the  sample  mean  as  a  test  statistic.   Letting 

if  Y.  >  Y 
J 

if   Y.  <  Y  , 
J 

for  j  =  l,2,...,n,  the  test  statistic,  Y,  can  be  expressed  as 


(4.1.1)  Y  =  E  cp.  . 

J-l  J 

Without  loss  of  generality,  we  can  assume  that  the  variables  in  the 

sample  have  been  standardized  to  have  zero  mean  and  unit  variance. 

Then,  using  the  representation  of  Y  given  by  (4.1.1)  and  assuming  that 

H  is  true,  we  can  find  the  mean  and  variance  of  Y.   For  the  mean, 
o  ' 

we  have 
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45 


[  A  "t  ] 


E(Y)  =  E 


n 
=   Z  E(cp.) 
3  =  1 


=     Z     P(<p.  =  1) 
3  =  1  3 


=      E  P(Y.  -  Y  >  0) 


n   i 
J  =  l 


n 
2  ' 


since  Y.  -  Y   is  symmetrically  distributed  about  zero.   Before 

calculating  the  variance  of  Y,  we  first  need  the  variance  of  cp  .  and 

covariance  of  c   and  CD  .   We  have 
J      k 

Var  (9.)  =  E(cp2)  -  [E(cp.)]2 
J       J        J 

-  E(CPj)  -  (|>2 

1 

~  4  ' 

and 

cov  (9j.<Pk)  =  E(cPjCpk)  -  E(tPj)  E«pk) 

=  P(Y.  -Y>0,  Y,  -Y>0)  -\   . 
3  k         4  . 

But  the  probability  above  is  identical  to  P_(p),  where  p  equals 

cov  (Y.-Y,  Y  -Y) 

corr  (Y.-Y,  Y,  -  Y) 
J      k 


v'Var  (Y .  -Y)  Var  (Y  -Y) 
J  k 
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The   covariance  of     Y.  -Y      and     Y,    -  Y      and   variance  of     Y.  -Y   are 

J  k  J 

cov    (Y.,Y,)-cov    (Y.,Y)-cov    (Y.Y,  )  +  cov    (Y,Y) 
J      k  j  k 

111 

=   0 +   - 

n       n       n 

1 

n 

and 

Var    (YJ  -   2   cov(Y.,Y)    +   Var    (Y) 

n       n 


-*-*• 


respectively.      Hence, 


P    = 


1 

n 


1-i 
n 


-1 
n  -  1 

Therefore, 


and 


cov   (cp .  ,9k)   =  Pa  (-  A.)  -  i 

1       .-If       1\ 

=    S—    Sln  I  -    T        . 

2tt  V     n-1/ 


Var    (Y)    =   Var    (   T,     cp .  J 


=      £     Var    (cp.)    +   2   E       I  cov    (cp    ,cp   ) 
j=l  J  j   <  k  J      k 


n       n(n-l) 

=    T    +    o SI 

4  2tt 


-1    /        1 
n        l"  5=1 
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The  standard  deviation  has  been  computed  for   n<  50   and  is  given 
in  Table  1. 

4.2   The  Null  Distribution 

Clearly  the  test  statistic,  Y,  has  a  discrete  distribution, 
taking  on  values   l,2,...,n-l  with  positive  probabilities.   The  proba- 
bility of  the  event,  [Y  =  r]  ,  r  =  l,2,...,n-l,  can  be  written  as 

Pn(Y=r)  =IP(CPil  =  1 ^r=1;C?ir+l  =  0 ^n  =  0)  ' 

where  the  summation  is  over  all  partitions  fi  i  ;  i   ,,....  i  } 

1      r   r+1      nJ 

of  the  set  [1,2 n] .   Then,  since  the  variables  q)  ,«p„ ,...  ,<P  are 

1      Z  n 

identically   distributed, 

P    (Y=r)    =    (*)  P(to    =1 cp    =1;   o         =  0,...,<P=0) 

n  \r/  l  r  r+1  n 

=    ()  P(Y„  -Y>0, Y     -Y>0;    Y         -Y<0, 

\r/  1  r  r+1 

...,Y     -Y<0)    . 
n 

After   the   transformation 


rr.    f  n   N5 


U.    =    (Y.  -Y)    (-=-) 
3  J  \n-lJ 


we  have 


P    (Y=  r)    a   (n)  p(u    >0, U    >0;    U      ,  <0,.  .  .  ,U    <0)  , 

n  Vr/  1  r  r+1  '    n 


where   the  normal   variables     D    ,U    ,,,,  ,U    ,    each  have  mean   zero   and 

1     <~i  II 

common  variance  and  correlation  given  by 


TABLE  1 
STANDARD  DEVIATION  OF  Y 
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n 


3 
4 

5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 


. 50000 

.59242 

.66760 

73389 

.79416 

. 84994 

.90214 

.95140 

.99818 

1.04283 

1.08562 

1.12678 

1.16646 

1.20484 

1.24202 

1.27811 

1.31320 

1.34738 

1.38071 

1.41325 

1.44505 

1.47617 

1.50664 

1.53651 


27 

1.56581 

28 

1.59456 

29 

1.62281 

30 

1.65058 

31 

1.67788 

32 

1.70475 

33 

1.73119 

34 

1.75724 

35 

1.78291 

36 

1.80822 

37 

1.83317 

38 

1.85779 

39 

1.88208 

40 

1.90607 

U 

1.92976 

42 

1.95316 

43 

1.97628 

44 

1.99914 

45 

2.02173 

46 

2.04408 

47 

2.06619 

48 

2.08806 

49 

2.19701 

50 

2.13112 

49 


Var  (U.)  =  -^—   Var  (Y.  -Y) 
J    n-1       j 


At    (1  "  -) 
n-1      n 


and 


=  1 


corr  (U.,U,  )  =  -^-  cov  (Y.  -Y,  Y,  -  Y) 
J   k    n-1       j      k 

=  Ar   ("  i) 

n-1    n 

1 
n-1 


Hence,  for  ra<n,  p  is  greater  than  -  — -  ,  so  that  the  m  variables, 

ra-1 

U, ,U_,...,U  ,  have  a  multivariate  normal  distribution.   Therefore, 
12      m 

P  (p  )  ,  k  <  m,  is  defined  and  can  be  computed,  using  the  method  discussed 

in  Chapter  3.   Furthermore,  since  P(Y=n)  =  0 ,  we  can  set  P  (p )  equal 

n 

to  zero  and  use  relation  (3.3.3)  to  find  P  (Y=r).   That  is, 

n 

(4.2.1)     P  (Y  =  r)  =  Of)     E   (-l)j  (nTr)  P   .(p)  , 

where       p  = and   P  (p )  -   0.   Using  equation  (4.2.1)  and 

n— l  n 

the  results  of  Chapter  3,  the  null  distribution  on  Y  was  obtained. 
The  results,  for  n  <  22 ,  are  given  in  Appendix  4. 


4. 3  Approximations  to  the  Null  Distribution 

David  (1962)  has  shown  that  the  asymptotic  distribution  of 

n  ( — )  is  normal  with  mean  zero  and  variance  —  -  —  .   Hence,  for 

\n  2/  4   2rr 

large  n,  we  should  be  able  to  approximate  the  distribution  on  Y,  using 
a  normal  distribution  function.   In  particular,  the  critical  values  of 
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Y  needed  to  form  the  rejection  region  can  be  determined  using  the 
approximation.   With  a  10%  level  of  significance  and  a  two-tailed 
alternative,  the  critical  values  are  the  solutions  to  the  equation 

1   n 

r+  2  "  2 

—————  =  ±  1.645  . 

VVar  (Y) 

(Notice  that  —  was  added  as  a  correction  for  continuity  factor.) 
As  an  example,  with  n  =  19,  we  have 

1    19 
P  +  2  "  T 
1.31320    =*1-M5 

or 

r  =  9  ±  2.1 
=  7,  11. 
As  a  check,  from  the  small  sample  distribution  on  Y  given  in  Appendix  4, 
we  have 

P19(Y-7)  =  P19(Yill)  =  -0595- 
Table  2  compares  the  small  sample  distribution  of  Y  with  the  normally- 
approximated  distribution  for   n=19.   (The  results  of  the  approxima- 
tion which  will  be  discussed  next  arelisted  in  the  third  column.) 
We  would  expect  the  approximated  results  to  increase  in  accuracy  as 
the  sample  size  increases. 

An  alternate  approach  to  the  problem  of  approximating  the  null 
distribution  on  Y  is  through  the  use  of  order  statistics.   Consider 

again  the  random  sample  of  standardized  variables,  Y„ ,Y„,...,Y  .   Then, 

12      n 

letting  Y ,  .  denote  the  vth  largest  order  statistic,  the  events  [Y<  r] 
and  [Y   .  <  Y]  are  equivalent  if  we  set  V  equal  to  n-r.   Letting  F  (x) 
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TABLE  2 
THE  CUMULATIVE  DISTRIBUTION  OF  Y,  n=  19 


Small 

Sample 

Distribution 

Approximations 

r 

* 
Normal 

Edgeworth'  sT 
Expansion 

1 

.0000 

.0000 

.0000 

2 

.0000 

.0000 

.0000 

3 

.0000 

.0000 

.0002 

4 

.0000 

.0000 

.0007 

5 

.0007 

.0011 

.0029 

6 

.0092 

.0112 

.0142 

7 

.0595 

.0639 

.0638 

8 

.2190 

.2231 

.2175 

9 

.5000 

.5000 

.5000 

10 

.7810 

.7769 

.7825 

11 

.9405 

.9361 

.9362 

12 

.9908 

.9889 

.9858 

13 

.9993 

.9989 

.9971 

14 

1.0000 

1.0000 

.9993 

15 

1.0000 

1.0000 

.9998 

16 

1.0000 

1.0000 

1.0000 

17 

1.0000 

1.0000 

1.0000 

18 

1.0000 

1.0000 

1.0000 

f 


Uses  first  two  moments  of  Y. 
Uses  first  four  moments  of  Y 


(n-r) 


-  Y. 
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represent  the  standardized  distribution  function  of  Y   -  Y,  we  can 
write 


(4.3.1)     P  (Y  <  r)  h  P(Y.  ,  -  Y  <  0) 
n  (v) 


=  P 


J  Y(v>-  Y  -  A 

I* 

2bV 


2*V  ' 


=  F. 


1 
2bV 


where   p.,  and  §   represent  the  mean  and  variance  of  Y   -  Y, 
1  V      2  V  ( v ) 

respectively.   An  approximation  to  the  distribution  F  (x)  can  be 
obtained  by  using  Edgeworth's  expansion  (see,  for  example,  Cramer, 
p.  229).   Letting  §    represent  the  kth  central  moment  of  Y   .-Y, 

K  V  \*) 

we  have,  using  the  first  four  moments, 


(4.3.2) 


Fv(x) 


=  G(x)  -  g(x) 


1    3^V 
6 


-3/2 
2S  V 


2 
(x  -1) 


where  G(x)  and  g(x)  are  the  standard  normal  distribution  and  density 

functions,  respectively.   Thus  letting   x  =  -  — s-  in  equation  (4.3.2), 

2*V 

we  can  approximate  P  (Y  ^  r) . 

n 

By  using  a  power  series  representation,  the  moments,   §  , 
k  -   2,3,4.  can  be  determined.   Saw  (1958)  has  shown  that  the  kth 


moment,   u  ,  of  the  Vth  order  statistic  can  be  expressed  as 

it  V 
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00 

(4.3.3)      ,  \k ,  =   E  H.(p..,0,k) 


j=0   j   V      (n+2)J 


=  E  H.(p,,,k)  n(j), 

where  d  =  and,  for  convenience,  we  have  replaced  H.Cp^.Ojk)  and 

*  v   n+1  J   v 

- — r  by  H.(p  ,k)  and  n(j),  respectively.   The  constants,  H  (p  ,k), 

(n+2)J      J   V  J 

are  tabul ated  for   k  =  1,2,3, 4 , j  <  5,  and  p^  =  . 50 ,  .  55 , 95  by 

Flora  (1965).   It  follows,  since  the  moments  ,  §.,  are  functions  of   k  , 

k  \>  k  v 

that 


k5V  =  .Z  Cj(Vk)  n(j)' 

where  the  constants   C.(p  ,k)  are  functions  of  H.(p  ,k). 

In  order  to  find  the  C.(pv  ,k),  we  must  first  express  the 

•J 

§   as  functions  of  the   u  .   Letting  9(t)  represent  the  character- 
istic function  of  V/v.-v~Yi  we  have 


Cp(t)  =  E 


[/«*(vf»] 


Since     Y   and  Y ,,   -Y   are   independent,    we  can  write 

(VJ) 


>(t)    =   E  le 


t(Y(v)-Y) 


r^tY] 


ELb' 


I6       J 


qt)2   M 

j 


2n  at) 

=  e  S       u 


j=o^v    J; 
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i=0  I]  (2n)        j=0  J  J' 


=      E        E      (-1) 


i   at)2^ 


j^   - 


£=0  j=0  i!j!(2n) 


After   the   change  of   variable,    k  =   2l  +   j ,    it   follows    that 


00     ,.itk 

■PCt)  -    =    ^L 

k=o      k- 


k! 


■k/2 

E        (-1)X  — 
Li=0  X!(k-2j£)!(2n) 


I   (k-2i)^V 


] 


Therefore, 


-  k  k/2  £  k' 

E(Y(V)-Y)       =    lA   +      E      (_1)         : 1    (k-2*A    ' 

(Vi>  kV        £=1  £!(k-2i):(2n)X    U  ^ 


and,    in  particular, 


E(Y(v)-  Y)  "  A   • 


E(Y(V)"  Y)2  =    2\   ~  £    • 


E(Y(V)  "7)3  =   3%   '  I  A    ■ 


E(Y(V)"Y)4     =  4^  *  S  2^v  +  4 

n 


Using  this  information,  the  central  moments,  ,  ^  ,  of  Y,   -  Y  can  be 
6  '  '  k  V      (v) 

expressed  as  functions  of   u   .   For  example,  with  k  =  2, 
2^  =   E(Y(V)-Y"  !'V2 

■  E(Y(vf Y)2  "  2E<Y(vfY)  A  +  1% 

1      2      2 

"2^   n    l^V  +  l^V 


2    1 

2  ►Si  "  l^V   n 
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Similarly , 

3§V  =  3'^V  "  Vv  1%  +  Vv  ' 

4§v  =  4^  +  |  1%  +  4  "  Vv  3^V  +  Vv  2% 

n 

6       _   4 
"  n  2^V  "  Vv  • 

Finally,  the  constants,  C.(p  ,k),  can  be  determined  as  functions  of 
H.(p  ,k)  by  substituting  the  power  series  representation  of  p,  and 
equating  coefficients  of  n(j).   With   k  =  2,  we  have 

F  2   I 

2"V   2^V  "  l^v  ~  n 

co 


r  T 

E  H  (p   2)  n(j)  -    E  H  (p   1)  n(j) 
j=0   J  U  j=0   J  J 

n+2  I       2  J 


1  "  n+2 


CO  00 

E  H  (p   2)  n(j)  -   E   H2(p  -1)  n(2j) 
j=0  J  j=0   J 


-  2  E  E  H  (p   1)  H  (p   1)  n(j+k)  -   E  2J  n(j  +  l) 


=  E  C  (p  ,2)  n(j) 
j=0  J 


Equating  coefficients  of  n(j),  we  find  that 


Co(V2)  =  W2)  "^V1*' 


C1(PV,2)  =  H1(pv,2)  -  Ho(pv,l)  Hl(pv,l)  -  1, 

C2(PV,2)  =  H2(pv,2)  -  H2(pv,l)  -  2Ho(pv,l)  H2(pv,l)  -  2, 
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C3(PV,2)  =  H3(pv(2)  -  2Ho(pv,l)  H3(pv,D 

-  2H1(pv,D  VPV'1}  -  4> 

C4(PV.2)  =  H4(pv,2)  -  l£(pv,l>  -  2Ho(Pv,D  H4Cpv,l) 

-  2H1(PV,1)  H3(PV>1)  -  8  , 

C5(PV,2)  =  H5(pv,2)  -  2Ho(pv,l)  H5(pv,l) 

-  2H1(PV,1)  H4(pvl)  -  2H2(pv>l)  H3(pv,l)  -  16. 

The  constants,  C.(p  ,3)  and   C.(p  ,4),  can  be  found  in  the  Same  tedious 
J   *  J 

manner.   The  values  of   C.(p  ,k),  k  =  2,3,4,  j<  5,  and 

p  =  .50,  .55,...,. 90,  are  tabulated  in  Tables  3,  4,  and  5. 

Using  these  tables,  we  can  approximate  the  moments  2§..>  3^v> 

and  §  ,  and  then  use  relation  (4.3.2),  with  x  = ,  in  order 

2^V 
to  obtain  an  approximate  distribution  on  Y.   As  an  example,  we  take 

n  =  19.   Then 


PV   n+1 

19-r 
"   20   * 

Tables  3,  4,  and  5  were  used  to  approximate  _§,_§,  and  §  , 

^r         2  v  3  v      4  v 

for  r  =  1,2,... ,9.   The  results  are  listed  in  Table  6.   The  values  for 

u  were  taken  from  tables  computed  by  Teichroew  (1956).   Teichroew' s 
tables  were  also  used  to  check  the  accuracy  of  the  series  approximation 
for   §  .   For  n  -   19 ,  the  approximation  was  accurate  to  five  decimal 
places. 
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These  moments  were  used  in  the  expression  for  F  (x)  given  by 
(4.3.2)   to  obtain  an  approximate  distribution  on  Y  for   n-19. 
The  results  are  given  in  Table  2,  page  51.   For   n =  19 ,  there  appears 
to  be  little  difference  in  accuracy  between  the  two  approximating 
methods — certainly  not  enough  to  justify  the  extra  labor  involved  in 
computations  for  the  latter  method.   However,  the  second  method  does 
work  well  and  is  at  least  of  theoretical  interest. 
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TABLE  6 


MOMENTS  OF  Y    -  Y,  n  =  19 


r 

* 
V 

** 

2^V 

§  t 
3^V 

4V 

1 

18 

1.37994 

. 11015 

.01897 

.04152 

2 

17 

1.C9945 

.07308 

.00868 

.01768 

3 

16 

.88586 

.05484 

. 00492 

.00975 

4 

15 

.70661 

.04416 

.00308 

. 00624 

5 

14 

.54771 

.03739 

. 00202 

.00442 

6 

13 

.40164 

.03298 

.00131 

. 00342 

7 

12 

.26374 

.03020 

.00079 

.00285 

8 

11 

. 13072 

.02866 

. 00038 

.00255 

9 

10 

. 00000 

.02816 

.00000 

.00246 

V  =  n  -  r. 


** 


L^  =  E(Y(v)-y). 


+  k^  =  E(Y(v)-Y-  i»Vk.   k=  2>3'4 


CHAPTER  5 


OTHER  METHODS  OF  EXPRESSING  P  (p ) 


5.1   Introduction 


Before  discovering  the  expression  for  P  (p  )  involving 

m 

Tcheby cheff-Hermite  and  Legendre  polynomials,  three  other  methods  of 

expressing  P  (p )  were  used  in  attempting  to  obtain  numerical  results. 

Each  of  the  first  two  methods,  outlined  in  Sections  5.2  and  5.3, 

expresses  P  (p )  as  an  infinite  series.   However,  in  each  case,  not 
m 

only  is  the  series  slow  to  converge,  but  no  workable  expression  can  be 

given  for  the  kth  term  of  the  series.   Therefore,  these  methods  are 

not  useful  in  obtaining  accurate  numerical  results.   In  Section  5.4,  we 

give  an  expression  for  P  (p  )  involving  the  moments  of  extreme  order 

m 

statistics.   However,  this  expression  can  be  used  only  for  limited 
values  of  m  and  rho. 


5.2  A  Power  Series  in  Rho 


Using  the  definition  of  P  (p )  given  in  Chapter  1,  we  can  write 

00        oo 

Pra(p)  =  f  ...  ,"  g(x  ,.  .  .  ,x  ;  p)  dx  ... dx  , 
m      J     ^     l      m       1     m 
o     o 

where  g(x  ,...,x  ;  p)  is  the  multivariate  normal  density  on  the  equi- 

correlated  variables  X  , .  .  .  ,X  .   Since,  when  rho  equals  zero,  the 

1      m 

density  function  simplifies  to 
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63 

i  ;  2 

m  -  -  Ex. 

"2    2  1=1  J 

g(x  ...  ,x  •  0)  =  (2tt)    e    J 
1      in 

_  co<x  <oo,  j   _  1,2,...  ,m, 
we  could  simplify  the  integrand  by  expanding  g(x  ,  ...,x  ;  p)  in  a 
Maclaurin  series  in  rho.   We  have 

(5.2.1)     g(Xl xm;  p)  =  =  6r(3g-)  lUj. V  p)|  n  ■ 

k=0      r  |p=0 

Unfortunately,  as  a  function  of  rho,  g(xH,...,x  ;  p)  is  quite  compli- 

1      m 

cated  and  it  is  not  feasible  to  take  derivatives  with  respect  to  rho. 
However,  the  following  identity  simplifies  the  problem  to.  some  extent: 

5^k 


(5.2.2)      (j-)  g(x1(...,xm;  p) 


p=0 

k 


i<J    i   3 
In  proving  identity  (5.2.2),  we  use  the  characteristic  function, 

cp(t,  ,...,t  ),  of  the  variables  X,  , .  .  .  ,X  .   By  definition 

1  m  1  m 

|    iE    t.x.j 
(5.2.3)  cp(t    ,...,t   )   =   E\e     J=1  / 


1)  •  •  • > ~    / 
m 


m 

00       iZt.x. 
J   J 


r      r      j=i  J  J 
=  j  •••  J    e  e<xi v  p)dxi- 


dx 
m 


—00  —00 


-i  t'vt 

e  E  —     — 


1      ™      ,2 


-s:    E     t.   -  p   E  E  t.t. 
.  2  J-l     J  i<J     x  3 


Differentiating  with  respect  to  rho,  we  have 
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5 
3p~  t 


00        _O0 


I 


m 
I    Z  t.x. 


_C0        —00 


g(x   .  ..  ,x  ;  p)dx    .dx 
1      m      1     m 


=  -  Z  Z  t  t .  e 


-  i'Z  t2  -  p  Z  Z  t.t 
2     J      ±<J   J  k 


Since,  for  -  — -  <  p  <  1,  g(x.,  ,.  .  .  ,x  ;  p)  is  a  continuous  function 
m-1  1      m 


5 


in  p ,  t—  g(x x  ;  p)  exists  and  is  integrable.   Therefore,  differ- 

cp    1     m 

entiation  inside  the  integral  is  permitted,  and,  at  the  point  p  =0, 
we  have 


(5.2.4) 


m 


!~6 

Urg(xi' 


,x  ;  p) 
ra 


I    Z  t.x. 
J  J 


1         j=l  3   3 

Je         dx^.. 


dx 


m 


p=0 


=  -  Z  Z  t.t .  e 


1  m   2 

-i  Z  t 

2  i 


i<j 


i  J 


Next,  we  consider  the  characteristic  function  of  X ,X  when  p  =  0. 

1      m 

From  (5.2.3)  we  have 


J 


m 
m  L    Z  t.x. 

J=l  J  J  ( 
e  g(x1,.. 


,x  ;  0)dx  . . . dx 
m      1     m 


=  e 


-1°   t2 
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Differentiation  with  respect  to  t  and  t  gives 

a      b 

m 


oo  L    E   t  X. 

1   6   5     j=l  J  J 

J  ST-'StT  e 

-°°   a   b 


g(x1, 


,x  ;  0)dx  . .  .dx 
m      1     m 


-1st2 

6   6     2  L     j 

3t-FtTe 

a   b 


so  that 


m 


00     °°        i  I  t,x, 

r         r  i=i J  J 

(5.2.5)         ...     (-xxj  e  J      g(x  ,...,x ;  0)dx  ...dx 
J      J      a  b  l      m      l     m 

_00  —00 


=  t  t,e 
a  b 


■i  E  t2 
2  >1  J 


Again,  differentiation  inside  the  integral  was  permitted,  since  the 

function  x  x,  g(x....,x  ;  0)  exists  and  is  integrable.   Since 
a  b    1      m 

Vb  gUl V  0)  ~   Fx-^rg(Xl'"-'Xm;  0)' 

a   b 


after  summing  both  sides  of  equation  (5.2.5)  over  values  of  a  and  b 
such  that  a  <  b,  we  have 


(5.2.6) 


•••  J  L 


m 


i.    E  t  .x. 
6   6,  1    j=i  J  J 


JJ-*r-xrs(xv--\;  0)J 

oo  ^a-^   a   b  ~* 


dx. . . . dx 


m 


=  E  E  t  t   e 
_   a  b 
a<b 


-I  "t2 
2  >i J 
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Finally,  addition  of  equations  (5.2.4)  and  (5.2.6)  gives 

00        CD 

/."•J    U-g(xi'""Vp)|    , 

-co    _co    r  'p=0 

"  *A  3ir*R:g(xi'---'xm;0)J 

a<b    z   b 

m 

■i   Z  t.x. 
J  J 


e  J  1      dx„ . . .  dx  =  ( 

1     m 


which    implies   that 


^g(xir...,xm;p)j        3"3T5Tg(Xi V0) 

'p=0      a  >•  b        a       d 


Since  derivatives  of  all  orders  exist  and  are  continuous,  the  preced- 
ing process  can  be  repeated  any  number  of  times.  .  Hence,  the  identity 
(5.2.2)  is  proven. 

Using  this  identity  in  the  Maclaurin  expansion  of  g(x,,...,x  ;p ) , 

1      m 

given  by  (5.2.1),  and  substituting  the  resulting  expression  into  P  (p ) , 

m 

we  have 

CO  00  00 

k                                      k 
(5.2.7)  Pm(p)    .  [..]    I    "-  (   SI     „!_    »    )     g(Xi V0)d^...d^ 

°       °  k=0  a     b        a       b 

GO  CO  CO 

klJ'-'J  {**   3ir'3T-)    e<x! V0)dxi"'dV 

k=0  °       °     a  <  b         a     T» 

The  utility  of  this  expression  for  P  (p )  depends  on  how  readily  each 

m  ' 

term  in  the  series  can  be  determined  and  on  how  quickly  the  series 
converges.   Differentiation  and  integration  are  no  problem  once  the 
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'  6   6  \ 

sum,  (E  E  r r — )  ,  has  been  expanded.   In  fact,  the  jth  derivative 

\  -s  v x  ^  x  / 
a<b   a   ta 

of  g(x  ,  ...,x  ;  0)  with  respect  to  x  is  given  by 
X      m  a 

(Sx")j  «<V-'V  0)  =  Wa  g(V)  (of)'  g(Xa> 
a  a 

=  (-1)J  H.(x  )  g(x   .  .  ,x  ;  0), 
J   a     1      m 

where  H.(x)  is  the  jth  order  Hermite  polynomial  in  x.   For  example, 


H  (x)  =  1,     H  (x)  =  x   -  1, 

O  £• 


H^x)  =  x,     H3(x)  =  x  -  3x. 


It  follows  that 


(5.2.8) 


'  (s)    g(x)  dx  = 


df_1  t  ? 
dx7    g(x) 


J 


=  (-1)J  *  H    (x)  g(x)J 


=  (-D 


For  even  values  of  j,  the  (j-l)st  order  Hermite  polynomial  vanishes 
at  the  point  x=0.   Therefore,  when  expanding  the  sum,  we  need  to 


consider  terms  that  involve  only  derivatives  of  odd  orders. 

Denoting  t —  by  0  ,  we  can  represent  the  kth  power  of  the 
a 
sum  by  the  multinomial  expansion 


(e  E  6  6)  =  )  k"   ,  (65)  X  (66)  2. . .  (6   6  )  N  , 

\  ^,   a  by    L  k, .  ■ ,  .  k„T.    12      13        m-1  m 


where  N  =  (I  and  where  the  summation  is  over  all  integer  values  of 

N 

k  ,ki;>,.  .  .  ,k  ,  satisfying   E  k.  =  k  and  k.  >  0,  j  =  1,2,...,N. 
J-        «  j=l  J  J 

Using  this  expansion  and  the  value  of  the  integral  given  by  (5.2.8), 

we  can  determine  the  first  few  terms  in  the  series  expression  (5.2.7) 

for  P  (p).   Denoting  the  (k+l)st  term  in  the  series  by  S  .  for  the 
m  •   k 

first  two  terms,  we  have 


Q  00  00 


O 


=  Por I  ■■■  !  g<v*-'V  0)  dxr--dxm 


-(? 


m 
2) 


and 


OS  00 

n 1    P  !»  r1  1  >  ki  ko  kvr 

h-  Ir  J -J  Zrrnrr  <W   <V3>   •••<8»-1V 

o   o    1     N 


g(x x  )dx   .  .dx 

1      m   1     m 


Since   E  k.  =  1,  there  are  N  =  (  J  terms  in  the  sum.   Therefore, 
j=l  J  W 

since  X1(...,X   are  identically  distributed,  it  follows  that 
1      m 
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00  cc 

5        5 


si  =  p  ©/•••/ SxT'*x:g(Ji---'V  °>dVdxi 


o   o    1   2 

2  ,^m-2 


o 

'  P  W  L"  ^-J   (2) 


V2;   — 


lV"   m(2) 

=  P 


where  m,  .  will  be  used  to  denote  the  permutation  of  m  elements  k 

at  a  time. 

N 
For  the  third  term,  S  ,   E  k.  =  2.   Since  a  value  of  k 

3=1      3  J 

equal  to  two  would  result  in  an  even  powered  derivative,  we  need 

to  consider  only  values  of  k.  equal  to  one.   Furthermore,  these  two 

"ones"  must  be  assigned  to  two  of  the  exponents  in  the  quantity 

kl      k2  kN 

(^i^c.)   ^1^0)   ...(6   6  )    so  that  the  resulting  product  contains 
12     13        m-1  m  b 

four  distinct  5's.   The  first  "one"  can  be  assigned  in  N  =  (_)  ways. 

k. 

Then,  there  remains  (    J  couples,  (*io^i4)   .   containing  6's  dis- 
tinct from  the  couple    (^i-.^.,)  »  determined  by  the  assignment  of 
the  first  "one."   Hence,  the  second  "one"  can  be  assigned  in  I    J 

ways.   Since  the  "ones"  are  not  distinguishable,  there  are  [?)[?    )  5- 

possibilities  for  selecting  the  two  non-zero  exponents.   It  follows 
that 
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p_  2!    /nN   /m-2\   1  p   *  fi   6   6   5 
2  "  2!  1!1!  \2)    '  \  2    J    *2  J  "  'J  Sx~* Fx~' 55T" W~ 


"X4 


g(x x  ;  0)dx   .  .dx 

1      m      In 


24   L"  v^J  ^ 


2  ....mm... 

P_  f  1N,    (4) 

2  \2/     2 


n 


Although  one  might  hope  for  a  general  expression  for  S  ,  such 
hopes  diminish  after  evaluating  the  next  two  terms  in  the  series. 
With   k  =  3,  we  must  consider  assigning  either  one  'three"  or  three 
'ones"  to  the  N  k.'s.   (The  choice  of  one  "two"  and  one  "one"  results 
in  an  even  powered  derivative. )   The  assignment  of  the  three  "ones" 
can  result  in  two  types  of  products  with  only  odd  exponents  on  the  6's. 

The  first  type,  6i16i26i36i4'  can  be  fo™ed  by  f^)  [2(m-2)  ]  (m-3)-^ 
different  assignments,  and  the  second  type,  6^  5^  b^   6^   6^   b±        can  be 

X    —    «J    4    O    D 

J  ^  fm\    'm-2Ym-4N  1    ^'^  „ 

formed  by  \^J       2  A  2  y  31"  dlfferent  assignments  of  the  three  ones  . 

Consequently, 

3  f       °°   co   3   3 

_    p   I  3!  /m\  f»    r>  6    6 

s3  =  3T\3T  [2)  J ■  •  -J  sr-^r  e(xr  •  •  •  -  V0)dxr  •  -dxm 

^  o         o         1         2 

aa  co         3 


+  irTnr(™)[2<"-2»<»-3>  £/••■/ 


6"       6        6        6 

6x„    5x„   6x„   5x„ 
o        o        1        2        3        4 


g(x    ,. . .  ,x    ;0)dx      .  .dx 
1  m  1  m 


1:1:1:  \2A  2  A  2  ;  s:  J' ••J  5"x" 

c 

....dx     1 
1  m  / 
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CO  CO 

6        6        6        6        6        6 


6x      6x„   0x„   6x„    6x^    Ox" 
ool23456 


g(x    , . . . ,x    ;0)dx 
1  m 


P3f(2)r    H2(0)-|2mm-2 


*1    *    L"   v^J    V"y 


+  m(4)   I 


r  V^ir   Ho(0h3/iNm-4 


m(6)  r  Ho(0h6m-8 


P_f  mm  r!m   !^n   m(6)i 

"   31     \2J     L    TT  2       +        3  J 


TT  TT 


Finally,  for  k=4,  values  of  the  k.'s  equal  to  "three"  and 


"one"  or  "two"  and  two  "ones"  can  both  be  assigned  to  form  a  product 

3   3 

having  the  form  ^i  6.  5.  5.  <   An  assignment  of  four  "ones  can 

X    m         O    4 
3 

result  in  either  a  ^i^i^i  ^i-Si  Si  product  or  a  6^6-^.5^ 
product.   It  follows  that 

„4  fT  /i.   /  \  I  r  Ho(0)-i2r  H  (0)-.2  ,  .m-4 

•  (T 
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+  Li:  1:1:1:  V2A  2  /I2  / 


or  Ho(0)i8  fi\m~ 


a- 


V^TT 


■]  (IT 


4     ,  „ .  m  r40m  .  . .  4m , „ . 

p_  /IN     1          (4)  (6) 

4;     \27     L        2  3 

tt  17 


It  is  unlikely,  from  the  expressions  given  for  S„  and  S  ,  that 
a  simple  general  expression  for  the  kth  term  exists.   Although  addi- 
tional terms  in  the  series   could  be  determined,  it  was  seen  from  numer- 
ical examples  that  the  series  converges  quite  slowly,  especially  for 
large  values  of  rho.   Some  of  these  results  are  given  in  Table  7. 


TABLE  7 


ERROR  INVOVED  IN  COMPUTING  P  (0 )  WHEN  THE  SERIES 

m  ' 

IN  RHO  IS  TRUNCATED  AFTER  FIVE  TERMS 


m=  5 

m=10 

P 

Absolute 
Error 

S4 

P 

Absolute 
Error 

S4 

-1/10 

.00001 

.00006 

-1/15 

.00006 

.00002 

-1/5 

.00008 

.00101 

-1/10 

.00029 

.00008 

1/5 

.00048 

.00101 

1/5 

.00318 

.00127 

1/10 

. 00003 

. 00006 

1/10 

. 00022 

.00008 

It  should  be  noted  that  this  method  of  computing  P  (0 )  is  valid  for 

m 

negative  values  of  rho.   In  fact,  if  the  results  in  Table  7  are  any 
indication  of  the  general  behavior  of  the  series,  we  would  expect  the 
fastest  convergence  for  p  <0. 
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5. 3  A  Series  Resulting  from  an  Inverse 
Taylor  Series  Expansion  of  g(u) 


For  oositive  rho,  consider  the  expansion  for  P    (p )  given 

r :  m 


in  Lemma  1: 


CO 

P   (p)  =  (m)  J     [1  -  G(9y)]r  [G(ey)]m"r  g(y)dy 
r :  m      vr  /  "  _ 


2         -1 
where   0  =  p(l-p)   .   Since  P   (p )  =  P     (p )  ,  we  can  write 

r : m       m-r : m 

P  (p  )  =  Pn   (p  )  ,  so  that 
m       0:  m 


P  (p)  =  /  G(0y)m  g(y)dy 
m      °  _ 


2 

After  the  change  of  variables,  u  =  8y  and  h  -    1/9  ,  we  have 

,  21 

(5.3.1)      p(p)=Ir  G(u)m  —  e     6  dy 

9  -=°       ^ 

=    (2tt)  2  h2  J     G(u)m  g(u)h-1  dG(u)  . 


By  expanding  the  density,  g(u) ,  in  an  inverse  Taylor's  series,  the 
integrand  will  contain  only  terms  involving  G(u) .  Expanding  about 
an  arbitrary  point  G(§) ,  we  have 


r«?  •*  *\             cr(^        V  CG(")  -  G(g)]J  (     d  V     I 
(5.3.2)      g(u)  =  I     j, (__,  g(u) 

j=0  ' 


=  J,  a (5)[G(u)  -  G(§)]J 


J 
j=0 
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The  function     a.  (u) ,    j  =  0 , 1, . . . ,      is  given  by 


1  (     d   \J 

<VU)  -  jy  VdGCuTy  g(u)- 


and  can  be  written 


jl    '        j!    VdG(u)/  L  dG(u)J 


1_  f     d     Y3"1  rdg(u)    /dG(u)~[ 
j.*    VdG(uW  L    du     /      du    J 


j:  VdG(u>y      u  ■ 


The  function 


[(dGTul)J  u]g(u)j  =  ~  e<u>J<J+i>:  Vi(u) 

is  tabulated  by  Saw  (1958)  for   j  =  0,1,..., 10. 

Substituting  the  expansion  (5.3.2)  of  g(u)  into  (5.3.1), 
we  have 

h-1 

BEJ 

h*  r   (  i  L 


1  ,  °°  .vh-1 


P  (p)  =  (2TT)  2   h2  f  (  E  a.(§)[G(u)  -  G(§)] 

o   j=0 


•  G(u)m  dG(u) 

h-1 

=  (2tt)  :   hS  J     z     0  (§;p)[G(u)  -G(§)]kG(u)mdG(u), 
o  k=0 

where  the  constant  3,  (§ ;p ) ,  k  =  0,1,...,  can  be  determined  from  the 
or  (§)'s.   Assuming  that  an  interchange  of  integration  and  summation 
is  permitted,  we  have 
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00 


Pi  f*  1cm 

P    (o)   =    (2tt)  h£      £     8,  (§;p)   J       [G(u)  -G(§>]    G(u)    dG(u) 

m  ,     „      k  J 


k=0 


Next,    defining     y^^'111)    by 

l 

yk(§,m)    =    \     [G(u)  -G(§)1      G(u)1"  dG(u)  , 
o 

we  can  write 

h-1 

00 

P   (p)   =    (2tt)    2     h^     S     6    (§;p)   y    (§,m)    . 
m  k=0     k  K 

A  recurrence  relation  f  or  y.  (§,m)  can  be  found  by  first  writing  the 

A 

integrand  as 

G(u)m[G(u)  -G(§)]k  =  G(u)m+1  [G(u)  -G(§)]k_1 

-  G(§)  G(u)m  [G(u)  -G(§)]k_1  , 


and  then  integrating  by  parts.   We  have 

1 


Yk(S,m)  :=  /  G(u)"*1  [G(u)  -G(§)]k_1  dG(u) 


-  G(§)  /  G(u)m  [G(u)  -G(§)]   X  dG(u) 


i  G(u)m+1  [G(u)  -G(§)]k 

~  o 


t 


1 

^~-  J      G(u)m  [G(u)  -G(?)]k  dG(u) 
o 


-  G(§)  Yk.!^-1") 


It      m_i_  "t 

=  k  [1  -  G(§)]   -  2±±  Yk(^"»)-G(§)  Yk_1<5,m) 
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It  follows  that 


(5.3.3)      Yk<5,m)  =  — ^  j_[l-G(§)]k  -  kG(§)  Y^^?^)]  • 


Since 

1 
Y0(f,m)  =  /  G(u)m  dG(u) 
o 

1 

"  m+1  ' 

the  Yk^'m)'  k  =  1,2,...,  can  be  obtained  easily  from  (5.3.3;  for 
given  values  of  m  and  G(§). 

The  selection  of  the  point  G(§)  should  be  made  so  that  the 
Yk(§,m)'s  are  small.   Since   (m+1)  G(u)    represents  the  density  on 
G(u) ,  we  can  write 

Yi/S.m)  =  -~  E[G(u)  -G(§)]k  . 
k        m+1 


Therefore,  by  letting 


G(§)  =  E[G(u)] 

=  /  (m+1)  G(u)m+1  dG(u) 


m+1 
=  m+2  ' 

Y1(5>m>  equals  zero  and  Y2^,m^  is  minimized.   However,  numerical 
work  has  shown  that,  although  Yk(5,m)   becomes  quite  small  as  k 
increases,  the  values  of  a (§) ,  j  =  0, 1, . . . ,   corresponding  to  the 
point  §  that  satisfies   G(§)  =  -~,  become  exceedingly  large. 
(For  example,  with   m=8,  §  =  1.28155157  and  a      (£)   =  -3,825,025.96.) 
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Therefore,  it  would  be  better  to  choose  the  point  G(§)  that  minimizes 
&.(§),  J  =  0,1,...,   especially  since  the  y,(§>m)'s  a**e  bounded  by 
(m+1)   .   Since  a\(§)  can  be  represented  by 

a   (§)  =  g(S) 
o 

^(5)  »  -  g 

:    i    J"2     i 

J      j!g(§)J   i=0  J>1 
where  the   a.  .   are  constants,  satisfying 

a.  .  =  0,         i  +  j  odd, 
a    <  0,         i  +  j  even, 

the  choice  §  =  0   clearly  minimizes  \a.  (§)|.   Denoting  a-.(O)  by 

j  j 

a. ,  we  have 
J 


ill 

(2tt)  2 

Of.  =  — - a  . 

J     3-  J,o 


so  that  a.   vanishes  for  odd  values  of  j.   Since  3  (0;p),  say  g  (p ) , 

u  K  K 

is  a  product  of  the  o-.'s  such  that  Zj.  =  k,  it  follows  that   3  (p ) 

3  i  k 

vanishes  for  odd  values  of  k,  and  that  the  series  P  (p )  contains  only 
even  terms. 

In  order  to  determine  the  values  for  the  o/.'s,  we  first  need 
the  coefficients  of  u  in  the  function  a.(u).   Denoting  this  coeffi- 
cient  by  cc    and  using  the  recurrence  relation  given  by  Saw  (1958), 
we  have 
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a         =  (2tt)  2 
0,0 

a.    .    =    .  ,2n, ;  [j-2)(j-3)a.  „   _ 
j,i   j(j-l)  J-2,1-2 

+  (2ij  -  5i  +  j-3)a.    . 

+  (i+l)(i+2)aj_2i_2],   i  <  j-2,  j  =  2,4,..., 

where  or.  .   equals  zero  if  either   i  >  j-1  or   i  <  0.   Using  this  rela- 
tion,  the  values  of  a,,    j  =  0,2, ...,22,  were  calculated.   They  appear 
in  Table  8. 


TABLE  8 


VALUES  OF  a . ,    j   =  0 , 2 , . . . , 22 


j 

a . 
J 

j 

a 
J 

j 

Of. 
J 

0 

.3989422803 

8 

-1.958451122 

16 

-105.6166131 

2 

-1.253314138 

10 

-4.703578753 

18 

-326.3330223 

4 

-.6562337483 

12 

-12.48581643 

20 

-1037.319292 

6 

-.9620889240 

14 

-35.44811307 

22 

-3373.253924 

Substituting  G(§)  =  G(0)  =  -  into  the  recurrence  relation 
(5.3.3)  for  Yi,(5,m)  and  denoting  y    (0,m)  by  y    (m)  ,  we  have 

,k 


,      %  1  f1^  k  /      xl 

Yk(m)    =   mTkTlL   V2J      "   2   Yfc-x^J 


Since  we  need     Yk(m)      only   for   even  values   of   k,    we  can   eliminate 
Yk_1(m)    in  the  above  relation.      It   follows   that 
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,  n        1  !~m    k(k-l)     m.%1   v   9  4 

Vm)  -  o^TITo^o  L^"1"- 4-Yk-2(m)J  •  k=  2'4---  • 

This  relation  can  be  easily  programmed  to  evaluate  v  (m)  ,  k  -   2,4 

for  a  given  value  of  m. 

Numerical  examples  were  used  to  investigate  how  quickly  the 

series, 

h-1 

l  °° 

(5.3.4)      P  (p)  =  (2tt)  2  h2   E  3,(p)  Yu(m). 

m  k=o  k    k 

converges.   As  examples,  we  include  the  numerical  computations  for 
m=10,  p=l/3,  1/4.   The  values  of  y.(rn),    3  (1/3)  ,  and  3  (1/4)  for 
k  s  0,2,..., 22   are  given  in  Table  9.   In  Table  10  is  listed 
h-1 

~T    i  n 

(2tt)    h2  E  3.  (p)  y.   (m) 
k=0  k     k 

for  n  =  0,2,..., 22  and  p  =  1/3,  1/4.   The  exact  value  of  P  (p )  is 

given  after  n =  »  .   As  can  be  seen,  absolute  errors  for  p  =  1/3  and 

p  =  1/4  are  .00258  and  .00122,  respectively.   In  addition  to  the  slow 

convergence  of  the  series  P  (p  ) ,  this  method  of  evaluating  P  (p )  has 

m  m 

two  other  disadvantages.   First,  the  expression  for   P  (p )  given  in 

m 

(5.3.4)  is  only  valid  for  p  =  - — — ,  h  =  1,2,...  .   Also,  the  constants 

h+1 

3  (p )  are  difficult  to  evaluate  for  small  values  of  rho. 

K 


TABLE   9 
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VALUES  OF  Yk<m).    B    (1/3) ,    0    (1/4) , 
k  =    0  ,  2  ,  .  .  .  ,  22  ,    m  =    10 


k 

Yk(10) 

8.  (1/3) 
k 

P.  (1/4) 

k 

0 

.0909090909 

.3989422803 

. 1591549430 

2 

.0163170163 

-1.253314138 

-1.000000000 

4 

.0032092907 

-.6562337483 

1.047197553 

6 

.0006629400 

-.9620889240 

.8772981706 

8 

.0001413557 

-1.958451122 

1.279624120 

10 

.0000308241 

-4.703578753 

2.418906535 

12 

.0000068352 

-12.48581643 

5.323901882 

14 

.0000015356 

-35.44811307 

12.95550088 

16 

.0000003486 

-105.6166131 

33.85865868 

18 

.0000000798 

-326.3330223 

93.33839437 

20 

.0000000184 

-1037.319292 

268. 1907545 

22 

.0000000043 

-3373.253924 

796.5361116 
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TABLE  10 


THE  FIRST  n  TERMS  IN  THE  SERIES  P  (p ) , 

m 


m=  10,  p  =  1/3,  1/4,  n  =  0,2, 


.22 


p  =  1/3 


t 


p  =1/4 


0 
2 

6 
8 
10 
12 
14 
16 
18 
20 
22 


.128564 
.056070 
.048604 
.046344 
.045362 
.044848 
.044546 
.044353 
.044222 
.044130 
.044062 
.044011 
.043753 


.157459 
-.020116 
.016459 
.022783 
.024757 
.025568 
.025964 
.026181 
.026309 
.026390 
.026444 
.026481 
.026603 


+  11" 

'Tabulated  entries  are  (2tt)2  •  2s   Z  3,  (l/3)v,  (10), 

k=0 


*  i   n 

Tabulated  entries  are  (2rr)  -3s   E  B,  (l/4)v,  (10) 

k=0  k      k 
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5. 4  Using  Moments  of  Extreme  Order  Statistics 

If  we  let  f(u;m)  represent  the  density  on  the  largest  order 
statistic  in  a  sample  of  size  m  from  a  normal  distribution,  then 


j.  ,        v    d  „ ,  .  m 
f(u;m)  =  —  G(u) 


=  mG(u)m    g(u),         _co<u<co. 


Using  the  representation  of  P  (p )  given  in  (5.3.1),  we  can  write 

m 

h-1 

_      (2tt)    h2  p°°   .h-1       ... 
Pm(p)  =  ~mTl J  g(u)    f^;1^1)  du- 


The  integral,  say  I,  above  can  be  simplified  by  successively  integrat- 
ing by  parts.   For  example,  after  integrating  by  parts  twice,  we  have 

.h-1  ._,  ,m+l 


I  =  /  g(u)n"  dG(u)' 


t   \n_1  ~,    Nm+1  I    r  ,u  ,>    ,    .h-1  „.  xm+l 
=  g(u)     G(u)      +  J  (h-l)u  g(u)     G(u) 


-oo   _eo 


du 


"  0+  S/ue(u)h-2  dG(u)' 


m+2 

—CO 


(h-l)(h-2)  f° r   2  ,x  ,      h-2  „,  vnw-2 


m 


>  U-ZJ  r  ,  2  _   .  .h-2     viih-2 
— ^ J  (u  -1)  g(u)     G(u)     du 


(h-1)  (h-2)  p°°  2  ,.   .  .h-3  .„,  ^m+3 

=  (m+2)(m+3)  J  (u  "1}  ^(u)     dG(u>     ■ 

_CO 

It  can  be  shown  by  induction  that  after  integrating  by  parts  k  times, 
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(h-l)(h-2).  ..  (h-k)   p  „  ,„,,  ^.^h-k-l. 


jy  J*  Hk(u)  gCu)11"11"  f(u;mfk+l)du, 


(m+2) (m+3).  .  .  (m+k+_. 


where  H  (u)   is  the  kth  order  Hermite  polynomial  in  u.   Therefore, 
k 


after  h-1  integrations,  we  have 

h-1 


2  ,  h „,   „  .  « 


(5.4.D       p  (p)  =  (2n),  i:;(h;;)(h;2);;-HJ  \  .ao  f(u;m+h)du 

m  (m+1) (m+2) ... (m+h)    J _  h-1 


—00 


h-1 
\  h  J        (m+h) 


where  U,    .is  the  largest  order  statistic  in  a  normal  sample  of 
(m+h) 

size  m+h. 

For  example,  with  p  =  1/3  and  p  =  1/4,  expression  (5.4.1) 
simplifies  to 

p   ^/3>   =  7 — iT7     ^    E  u 

m  (m+1) (m+2) 

(m+2) 
and 


p  u/4>  -         4  ^  n  T  e 

mv  (m+1) (m+2) (m+3)   L  „ 

(m+3) 


,2  - 1] , 


respectively.   Using  the  table  of  moments  of  extreme  order  statistics 

computed  by  Ruben  (1954),  expression  (5.4.1)  can  be  used  to  determine 

P  (p )  for  values  of  m  and  p,  satisfying  p  2:  1/12   and  m  +  p   <51, 
where  p    is  a  positive  integer. 
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APPENDIX   1 

Jk<p) 

k 

1/    2 

p 
1/    2 

1/  4 

1/    5 

0 

1.0CC0C 

l.COCCC 

l.CCCOC 

l.OCQCO 

2 

C.OOOOC 

-C. 17548 

-0.25871 

-C. 30772 

4 

O.OOCGC 

-C.C0677 

0.C2771 

G.C6523 

6 

o.ccccc 

-C.C0168 

C.CC243 

-0.C0299 

9 

C.QCCOC 

-C.CCC52 

C.CCC45 

-C.C0041 

10 

o.ccooc 

-C.C0C21 

0.CCC12 

-0.00008 

12 

c.oococ 

-C.CCC1C 

C.CCC04 

-0.COCO2 

14 

c.ccccc 

-C.CCCC6 

C.CCCC2 

-C.CCCCl 

16 

c.ncccc 

-C.CCCC3 

C.CCC01 

-C.COOCC 

1 8 

C.OCOGC 

-C.CCCC2 

O.CCCCC 

-C.COGCO 

2C 

o.ococc 

-C.00CC2 

o.ccccc 

-C.COOOO 

22 

o.occcc 

-C.CCCC1 

c.ccccc 

-O.COCCO 

1/6  1/7  1/8  1/    9 


c 

l.OCCGC 

l.CCCOC 

l.CCCOC 

i. COCCC 

2 

-0.34C1C 

-C. 36211 

-C.38C32 

-C. 39368 

<t 

0.19766 

C.  12483 

C. 14757 

C.16&73 

6 

-C.D1385 

-C.C2698 

-C.C4C7C 

-0.05415 

8 

0.OC002 

C. 00241 

C.CC647 

C.C1172 

LC 

C.0CCC4 

C.CCCC8 

-C.CCC32 

-C.C0134 

12 

C.CCC01 

C.CCCCC 

-C.CCCC2 

C.CC002 

14 

C.OCCOC 

-C.CCCCC 

-C.CCCCC 

0.00000 

16 

C.OCOGC 

-C  .COCCC 

-O.CCCCC 

C.OOGCO 

18 

O.OCCCC 

-C.CCCCC 

-O.CCCCC 

O.COCCO 

20 

O.OCOOC 

-C.CCCCO 

-C.CCCCC 

C.COOOO 

22 

O.OCCCC 

-C.OCCCC 

-C.CCCCC 

0.00000 

1/10  1/11  1/12  1/13 


0 

l.COCCC 

l.CCCCC 

l.CCCCC 

l.COOCC 

2 

-C.4C425 

-C.413C7 

-C.42C33 

-0.42647 

4 

0.18303 

0.197C4 

C.2C919 

C. 21981 

6 

-C. 06696 

-C. 07895 

-C.C9CC9 

-C. 10040 

a 

CO  1776 

C. 02425 

C.C3C96 

C.C3771 

10 

-0.0C3OC 

-C.CC524 

-C.CC796 

-0.01107 

12 

0.OCO24 

CCCC69 

C.CC140 

0.00239 

14 

C.CCCCC 

-C.CCCC3 

-C.CCC14 

-C.00G34 

16 

-C.OCCOC 

-C.CCCCC 

O.CCCOO 

C.C0C02 

18 

-C.OOOOC 

-C.CCCCO 

C.CCCOO 

0.00000 

2: 

-O.OCOOC 

C.CCCCO 

O.CCCOC 

-O.CGGCO 

22 

-C.OCCCC 

-C.CCCCC 

C.CCCCC 

-0. COCCC 

85 
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k  1/14  1/15  1/16  1/17 


0 

l.OOGOC 

l.CCCCC 

l.CCCCC 

1.00000 

2 

-0.43173 

-C. 43629 

-G.44C23 

-C. 44350 

4 

C. 2 29 17 

C. 23749 

C. 24491 

C.  25159 

6 

-0.1C992 

-C. 11872 

-0.12686 

-C. 13439 

3 

C.0444C 

C.05C96 

C.Q5733 

C.C6348 

1C 

-C. 01447 

-C.C18C9 

-C-C2185 

-C. 02570 

12 

0.00362 

C.C0512 

0-C0682 

C 00870 

14 

-C.CGC66 

-C.Q011C 

-C.CC168 

-C.C0239 

16 

O.0CCO7 

C.C0C17 

C.CCC31 

C.C0051 

18 

-cooooc 

-CCCCC1 

-C-CCCC4 

-0.C0008 

20 

-o.ooooc 

c.coccc 

C-CCCOO 

C.C0001 

22 

-c.occcc 

c.ccccc 

-c.cccoc 

-o.cooco 

1/18  1/19  1/2C  1/21 


0 

l.ccooc 

l.OCCCC 

l-CCCOO 

l.OOCCO 

2 

-0.44692 

-C. 44972 

-0-45223 

-C. 45451 

4 

0.25761 

C.263C8 

0. 26  80  7 

0.27263 

6 

-0.14136 

-C 14737 

-C-  15390 

-C. 15953 

8 

0.06941 

C. 07510 

C.G8C55 

C.0S577 

IC 

-C. 0296  1 

-C.C3352 

-0-.C3742 

-C. 04129 

12 

0.01C72 

C. 01289 

C-C1516 

C. 01751 

1^ 

-C.:C323 

-CC0419 

-C- CO 52 6 

-C. 00644 

lc 

0.OCC78 

C. 00112 

C-CC154 

C.002C2 

13 

-o.rcni? 

-C.CCC24 

-C.CCC37 

-C.C0053 

23 

0.CCOU2 

C.CCCC4 

0.CCCC7 

C.C0G11 

22 

-o.ocooc 

-C.COCCC 

-C-CCC01 

-C. 00002 

1/22  1/23  1/24  1/25 


n 

l.CCCCC 

l.CCCCC 

L-CCCOC 

1.00000 

2 

-0.4565? 

-C.  45847 

-C.46C2C 

-C. 46179 

4 

0.27682 

C  .28C68 

C.  28426 

C. 28757 

6 

-0.16479 

-C. 16971 

-C. 17433 

-0.17866 

3 

0.09076 

C.  09 5 53 

C-1CCC9 

C. 10445 

10 

-C.0451C 

-C.C4685 

-C.C5253 

-0.05613 

12 

0.0199  1 

C  .  C  2  2  3  7 

C-C2465 

C.C2735 

14 

-C.0C771 

-C.CC9C6 

-C.C1C48 

-0.01196 

16 

0.00258 

C.CC22C 

C.CC39C 

0.C0465 

18 

-O.OC073 

-CC0C98 

-C.CC126 

-C.G0159 

2  0 

0.00017 

C.CCC25 

O.CCC35 

C.C3047 

11 

-0.CC003 

-C.0CCC5 

-c.cccce 

-0.CCG12 
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-1/    3  -1/    4  -1/    5  -1/    6 


0 

l.OCCOC 

l.CCCCC 

l.CCCCC 

1.C00CC 

2 

-0.6245;. 

-C. 74129 

-C. 69226 

-C.fc599C 

A 

2.44596 

1.51C71 

1.15617 

C. 96395 

6 

-4.C4666 

-2.53547 

-1/    7  -1/    £  -1/    9  -1/10 


0 

l.OCCOC 

l.CCCCC 

l.CCCGC 

l.COCCO 

2 

-0.63689 

-C. 61968 

-0.6C632 

-C. 59565 

4 

0.35336 

C.775C5 

C. 71856 

0.67592 

6 

-1.86571 

-1.49215 

-1.2586C 

-1.09863 

S 

7.0626C 

4.46698 

3.2C847 

2.48636 

10 

-12.72389 

-3.09368 

-1/11  -1/12  -1/13  -1/14 


0 

l.OCCOC 

l.CCCCC 

l.CCCGC 

l.COCCO 

2 

-0.58692 

-C. 57967 

-0.57353 

-0.56827 

4 

C. 64262 

C. 61592 

0. 594C3 

C. 57576 

c 

-0.98317 

-0.89624 

-C.E2862 

-C. 77462 

5 

2.0272C 

1.71284 

1.46657 

1.32C03 

1C 

-5.7169C 

-4.32794 

-2.43337 

-2.83220 

12 

22.39377 

14.93534 

10.42892 

7.75038 

14 

-43.61756 

-27.917C0 

-1/15  -1/16  -1/17  -1/18 


0 

l.OCCOC 

l.OCCCC 

l.OCCOC 

l.COCCO 

2 

-0.56371 

-C. 55972 

-0.5562C 

-C. 55303 

4 

C. 56029 

C.547C3 

0.53553 

C. 52546 

6 

-0.72058 

-C.694CC 

-C. 66318 

-G. 63687 

8 

1.  18986 

1.08676 

1.CC33C 

0.93454 

10 

-2.39876 

-2.07684 

-1.83032 

-1.63672 

12 

6.03114 

4.86226 

4.C3107 

3.41332 

14 

-19. 31984 

-14.15336 

-10.82916 

-3.57463 

16 

82.16048 

52.68241 

36.20290 

26.21621 

18 

-155.97371 

-ICC. 15185 
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k  -1/19  -1/2C  -1/21  -1/22 


c 

1. 00000 

1.CCCC0 

1.C0C0C 

1.00000 

2 

-0.55028 

-C. 54777 

-C. 54549 

-C. 54342 

4 

.      0.51658 

0.50868 

0.5C162 

C. 49527 

6 

-C. 61415 

-C. 59426 

-C. 57696 

-0.56155 

8 

0.87702 

C. 82827 

0.78648 

0.7503C 

10 

-1.48144 

-1.25464 

-1.24948 

-1.16111 

12 

2.95302 

2.59C89 

2.30307 

2.C7017 

14 

-6.97997 

-5.31267 

-4.92343 

-4.255C0 

16 

19.77901 

15.42370 

12.3584C 

10.12926 

18 

-68.43842 

-49.08520 

-36.59633 

-28.16580 

20 

297.927CS 

191.51C78 

130.2634C 

92.67789 

22 

-571.91135 

-367.95599 

-1/23  -1/24  -1/25  -1/26 


0 

l.OCCOC 

l.GCCCC 

1.CCC0C 

1.C02C0 

2 

-0.54152 

-c.53<;ec 

-C. 53821 

-C. 53674 

4 

0.48951 

C. 48429 

0.47952 

C. 47514 

6 

-C. 54762 

-C.  53550 

-C. 52439 

-C. 51432 

8 

0.7187C 

C.69C69 

0.66622 

0.64422 

10 

-1.0E596 

-1.C2140 

-C. 96543 

-0.91650 

12 

1.87876 

1.71931 

1.58483 

1.47034 

14 

-3.72057 

-2.292CC 

-2.94291 

-2.65465 

16 

8.^6292 

7.18757 

6.19143 

5.39946 

13 

-22.2593C 

-17.98S9C 

-14.82C64 

-12.41346 

20 

63.4C097 

52.04764 

40.63786 

32.43603 

22 

-249.31169 

-176.16453 

-128.6  642  5 

-97.09741 

24 

1102.36011 

709.76567 

479.32931 

336.70025 

26 

-2132.Z2218 

-1373.66691 

APPENDIX  2 


m 


C 
2 
4 


C.  2333333  3 
0.  16666  66  7 


C.2  5C0CCCC 
0.25CCOCCC 


0.2CCCCC0C 
0.23571429 
0.28571429 


0.16666667 
0.29761905 
0.29761905 


8 


4 

6 
8 


u 
0 
C 


14285714 
29761905 
05844156 


0.CC10822: 


C.1250CCC0 
0.29166667 
0.07954545 
0.0C378788 


0.11111111 
0.28282628 
0.097  90210 
0.00808081 

0.CC308C81 


C.1C000CO0 
0.27272727 
0.11328671 
0.01363636 
C. 01363636 


10 


11 


12 


13 


0 
2 

4 

6 

8 

10 


O.C909C909 
0.26223776 
0.1258  7413 

0.0  20 j  5  34  £ 
O.C0G92C13 
G.CC00054  1 


0.08233333 
C.2518315C 
0. 13598901 
C. 02696078 
C.0C18556C 
0.CCC02977 


0.07692308 
0.24175824 
0.14396836 
C.C3405573 
O.C031S103 
O.CCC09319 
0.CCCC9319 


C. 07142857 
0.23214286 
C.15021CG8 
0.04111013 
C.0C488722 
0.00021874 
C.0CC21874 


.2 
4 

6 
8 
10 
12 
14 
16 


0. 

c. 

0. 
0. 
0. 

o  - 


14 

06666667 
22303922 
15495356 
04796182 
0C694127 
0C042872 
CCCCC872 
CC0OCO0  2 


15 

0.0625CCC0 
C.21446C78 
0.15847523 
0.054502C6 
0.00929634 
0.0CC742C2 
0.0CC02337 
0.CCCC0C19 


16 

0.C5S82353 
0.2C639835 
0.16O99C71 
0.C6C66317 
0.01189931 
0.CCH7258 
0.CCOC5157 
0.COCCCC77 
O.CCCCCC77 


17 

C. 05555556 
C. 19883041 
G.  16267943 
C. 06640778 
C. 01469616 
0.00172896 
0.0G009937 
C.0C00O232 
C.CC000232 
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m 
k  13  19  2C  21 

0       O.C5263158     C.05CC0CCC     0.04761905  C. 04545455 

2       0.19172932     C. 18506494     0.1788C670  C.  17292490 

4       Q. 16363936     C 16414455     C. 16414455  0.16377318 

6       0.07172C41     0.0766CC79     0.081C5904  C. 08511199 

8       0.01762539     0.02067C05     0.C2375868  0.02686559 

1C       0.C0241458     O.CC322t39     G.C0416566  C.0C521391 

12       O.CC017309     C.CCC27695     0.CCL42265  C. 00060912 

14       O.CC00G570     C.CCCC1211     0.0GC023C7  C.0CC04038 

16      0.CC00CC07     C.CCCC0C22     0.G0CC0C59  C.0C00C137 

18       0.CC00CC0C     O.CCCCCCCC     0.C0CCCC01  0.OCOOO0O2 

20  0.CCCCCC01  C.0C0C0C02 


22             23             24  25 

C  0.C4347326  C. 04166667  0.C4CCCC0C  C. 03846154 

2  0.16739130  C. 16217949  0.15726496  C. 15262515 
4  0.16309922  0.16217949  0.161C6101  0.15978275 
6  0.08378061  0.09208812  0.C95C5871  C.0977166C 

3  0.02996C7C  0.02301919  0.03602C93  0.03894997 
10  J. 00637367  CCC763419  0.CC897407  C. 01038666 
12  0.C003423Z  C. 00112522  0.CC145975  C. 00164685 
14  0.CC0r'6602  C.0CG10212  0.CCC15C82  C.CC021423 
16  0.CC0OC284  0.0CC00536  0.CCCCC942  0.0C001557 
lo  0.C00)O006  C.CCCC0C15  0.CCCCCC33  C.CC0C0066 
20  C.CC00U30C  O.OCCCOCCO  0.C0C0CC01  C.0CCC0G01 
22  C.GCO-CCOC  C.CtCCCCCO  O.CCCCGCOO  C.OCOCCCOO 
24  O.CCCCCCCC  c.ccccocco 


m 


APPENDIX  3 

P    (p) 

m 

1/2  1/3  1/4  1/5 


1 

0.33333 

C.304C9 

C.29C22 

C.282C5 

— 1 

0.250CC 

C.2C613 

C.  18532 

C. 17307 

4 

C2000C 

C.14S74 

0.12648 

C. 11301 

? 

0.16667 

C. 11413 

C.09C66 

C. 07741 

6 

0.14286 

C.C9C12 

0.06748 

0.05508 

7 

0.125CC 

C. 07311 

0.C5176 

0.04043 

8 

n.iim 

C.06C61 

0.C4C67 

C. 03044 

9 

0.1CC0C 

C.C5113 

C.C3262 

0.02343 

10 

C. 09091 

C. 04375 

0.0266C 

0.01836 

11 

0.08333 

C.C3790 

0.02-202 

C. 01463 

12 

C. 07692 

C. 03318 

C.C1845 

0.01182 

13 

0,0714  3 

C.02S3C 

C. 01564 

0.C0967 

14 

0.06667 

C.C26C8 

0.01338 

C. 00799 

15 

0.0625C 

C. 02338 

0.C1155 

0.C0668 

16 

C. 05882 

C.C21C8 

C.C1C04 

0.00563 

17 

0.05556 

C. 01912 

C.CG879 

0. 004 7 8 

13 

0.05263 

C. 01742 

C.CC775 

0.00^09 

19 

0.0  500C 

C.  01594 

C.CC687 

0.C0352 

20 

C. 04762 

C. 01465 

C.CCfcl2 

C.C0305 

21 

0.D4545 

C. 01351 

C.CC548 

0.00266 

22 

0.04343 

C. 01251 

0.C0492 

0.00233 

1/6  1/7  1/8  1/9 


2 

0.  27665 

C.272ei 

3.26595 

0.26772 

3 

0.16493 

0.15922 

C.15492 

0.15158 

4 

0.1G422 

C.09EC4 

C.C9345 

0.08990 

5 

C.06R92 

C.062C6 

C.C5375 

G. 05546 

6 

C. 04733 

C. 04205 

C.C3325 

C. 03538 

7 

G. 03352 

C.02892 

C. 02566 

0.02323 

Z 

0.02*37 

C.C2C42 

0.01766 

0.01565 

9 

C.  01612 

C  .  0 1 4  7  4 

0.C1244 

0.C1079 

10 

0.01374 

C.C1C86 

CCC894 

C.C0758 

11 

CO  105  9 

c.ooei4 

C.CC654 

0.00543 

12 

O.0O925 

C.0C62C 

CCC486 

C.C0395 

13 

C.0C653 

C.CC479 

C.C0366 

C.C0291 

14 

0.00526 

C.CC274 

C.C028C 

C.C0218 

15 

C. 00428 

C. 00296 

C. 00216 

C. 00165 

16 

0.0C351 

C.C0237 

C. 00169 

C. 00126 

17 

0.0C29C 

C.C0191 

C.C0133 

C. 00098 

18 

C.  00242 

C  .00155 

C.CC1C6 

0.COG76 

19 

C002C3 

C.0C127 

C.CC085 

0.C0C60 

20 

0.00172 

C.C01C5 

0.CCC69 

C.C0C48 

21 

0.00146 

C.CCC38 

C  .  C  C  C5  6 

C.C0C38 

22 

O.0C125 

C.C0C73 

C.CCC46 

0.CG031 

91 


92 


m 
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1/10  1/11  1/12  1/13 


7 

•— 

C. 26594 

C. 26449 

0.26328 

0.26225 

2 

0,14891 

C.14673 

0.  14492 

C.  14336 

4 

C. 08709 

C .08480 

C. 08289 

C. 08129 

5 

0.0  528  6 

C.05C77 

C.C49C4 

C. 04759 

0 

C  .  0  3  3 1 4 

C.03135 

0.02989 

0.02867 

7 

0.02137 

C.  01990 

C. 01870 

0.01772 

8 

0  .  0  1 A 1 3 

C. 01294 

C. 01198 

C. 01120 

9 

0.00955 

C.C086C 

0.CC784 

0.00724 

1C 

0. 00659 

C.0C583 

C.0C524 

0.00476 

11 

0.00462 

C.C04C2 

C.C0356 

C. 00319 

12 

0.0C33C 

C. 00282 

C.C0245 

C.C0217 

13 

0.00239 

C0C2C1 

0.CC172 

0.00150 

14 

COO  175 

C. 00145 

C.CC122 

C.C01C5 

15 

C.0C13C 

C.001C6 

C.CCC83 

C. 00074 

16 

0.0009E 

0.CCC7  8 

C.CC064 

0.00053 

17 

O.OC074 

C.C0C58 

O.CCC*7 

0.00039 

13 

0.0C057 

C.CCC44 

C.CCC35 

0.C0028 

19 

0.0C04  4 

C.C0C34 

C.CCC26 

G.C0021 

20 

0.00034 

C.CCC26 

C.CCC2C 

C.C0016 

21 

0.00027 

C.C0C2C 

C.C0C15 

C.C0012 

22 

0.0CC22 

C.CCC16 

C.CCC12 

C.00CC9 

1/14  1/15  1/16  1/17 


2 

C. 26138 

C.26C62 

C. 25995 

0.25937 

3 

0.14207 

C.14C93 

0.13993 

C. 13*05 

4 

C. 07992 

C.C7874 

C. 07771 

0.07680 

5 

0.04636 

C.0453C 

C. 04438 

C.04357 

6 

0.02764 

C. 02676 

C.C26CC 

0.02533 

7 

C. 01669 

C. 01619 

C.C1559 

C.C15C6 

3 

0.01C56 

C.C1CC1 

C.CC954 

0.00914 

9 

0.00673 

C.CC631 

0. 00596 

0.00566 

10 

G. 0043 8 

C.C04C6 

0.00379 

0.00356 

11 

0.00269 

C.0C265 

C.C0245 

C. 00228 

12 

0  . 0  C  1 9  4 

C. 00176 

C. CO  161 

C. 00148 

13 

0.00132 

C. 00118 

0.C0107 

C. 00098 

14 

0.0C091 

c.cocei 

C.CCC72 

0.00065 

15 

0.0C064 

C.C0056 

C.CCC49 

C.C0C44 

16 

0.0004  5 

C.C0C39 

C.CCC34 

C.C0030 

17 

C.0C032 

C.C0C28 

C.Z0C24 

C. 00021 

18 

0.0C023 

C.0CC2C 

C.CCC17 

0.00015 

19 

0.2CC17 

C.0CC14 

0.CCC12 

C.C0C1O 

20 

0.00013 

C.C0C10 

C.CCC09 

C.CC007 

21 

0.00009 

C.C0CC8 

CO  CCO  6 

0.00005 

22 

C.0C0C7 

C.C0CC6 

C.CCC05 

C.OCC04 
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m 
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1/18  1/ 19  1/2C  1/21 


2 

0.  25885 

C.25d38 

C. 25796 

C. 25758 

3 

0.12827 

C. 13757 

0.13694 

C.  13637 

t\ 

C. 07599 

C. 07527 

C.C7462 

0.C7403 

5 

C. 04285 

C.  04222 

C.C4165 

C.  04113 

6 

0.02475 

C. 02423 

C. 02376 

C. 02335 

7 

0.0146C 

C. 01420 

0.C1334 

C. 01352 

8 

C.0C879 

C.  00849 

'  C.CC821 

C.C0797 

9 

C.0C539 

C  .CC516 

C.CC49  5 

C..C0478 

10 

0  .  0  C  3  3  6 

C.CC22C 

C.CC3C5 

C. 00292 

11 

O.OC212 

C. 00201 

C.C019C 

C.  00180 

12 

0.00137 

C. 00128 

0.C012C 

C.C0113 

12 

C.0CC9C 

C  .CCC83 

C.CCC77 

C.C0G72 

14 

0.GCC59 

C.CCC54 

C.CCC5C 

0.CC046 

15 

0.0C04C 

C.C0C36 

C.CCC33 

0.C003C 

16 

0.C0C27 

C.C0C24 

C.CCC22 

C.CC02C 

17 

C.00C18 

C.CCC16 

C.CCC15 

G. CO 01 3 

18 

0.0CC12 

C.CCC11 

c.ccci: 

C.C0CC9 

19 

C.0CCC9 

c.cccca 

0.CCCC7 

C.C0G06 

2G 

C.0C0C6 

C.CCCC5 

0.CCC05 

0.C0G04 

21 

C.0CCC4 

C.CCCC4 

C.CCC03 

C.C00C3 

22 

C.0CCC3 

C.CCCC3 

C.CCCC2 

C.C0CC2 

1/22  1/23  1/24  1/25 


2 

0.25724 

C. 25692 

0.25663 

0.25637 

3 

0.13586 

C. 13538 

C. 13495 

C. 13455 

4 

C.0735C 

0.072C2 

C.C725S 

C. 07217 

5 

0.04C67 

C.C4C24 

C.C3985 

C. 03950 

6 

0.02297 

C. 02263 

C.C2232 

0.02203 

7 

0.01322 

C.  01296 

0.01273 

0.01251 

Q 

C.CC775 

C.C0756 

C.CC738 

0.00722 

9 

C.0C462 

C.C0448 

C.CC435 

C. 00423 

10 

C.0C28C 

C.CC27C 

C.CC26C 

C.C0252 

11 

0.0C172 

C. 00165 

C.CC158 

0.C0152 

12 

0.00107 

C.001C2 

C.CCC97 

0.C0093 

13 

0.00068 

C.CCC64 

C.CCC61 

0.C0053 

14 

0.CC042 

C.CCC41 

0.CCC38 

C.0OG36 

15 

0.00028 

C.CCC26 

C.CCC24 

C.CC023 

16 

0.00018 

C.C0C17 

C.CCC16 

C. 00015 

17 

CO  00  12 

C.CCC11 

C.CCC1C 

C.C00C9 

18 

O.0C0G5 

C.CCCC7 

C.CCCC7 

G.00G06 

IS 

0.0CCC5 

C.CC105 

C.CCCC4 

C.C0G04 

20 

0.OCCC4 

C.CCCC3 

0.CCC03 

G.C0G03 

21 

0.0CCC2 

C.C0CC2 

C.GCC02 

C.C00C2 

22 

C.0CCC2 

C.CCCC2 

C.CCC01 

G. 00001 
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APPENDIX  3  (Cortt 

int 

led) 

n 

-1/  2 

-1/  ? 

P 

-1/  4 

-1/  5 

2 
3 
4 
5 

0.16fc67 

C.  19 591 
C.G42E7 

C.2C978 
0.C6468 
C.CC978 

C.  21795 
0.07693 
0.01872 
C. 00192 

-1/    6  -1/    7  -1/    8  -1/    9 


2 

0.22335 

C.  22719 

C.23CC5 

0.23228 

3 

C.D85C2 

C.C9C78 

0.09508 

0.09842 

4 

0.~>253C 

C.C3C22 

0.C3402 

0.C3703 

5 

C  .  0  0  4  8 1 

C.C0759 

C.CC992 

C.C1168 

6 

C.0CC3M 

C.CC116 

C.CC211 

C. 00304 

7 

C.CCCC6 

C.CCC26 

C.C0C55 

8 

C.CCC01 

0.C0C05 

9 

O.OOOCC 

-1/10  -1/11  -1/12  -1/13 


2 

0.22406 

C. 23551 

C. 23672 

0.23775 

3 

0.10109 

C.1C227 

C.1C5C8 

0.10662 

A 

0.03947 

C. 04149 

0.C4318 

0.04462 

5 

C. 01352 

C. 01494 

0  .  C 1 6 1 4 

C 01719 

6 

0.00389 

C. 00467 

C.C0536 

G.C0b98 

7 

C. OOC37 

C.CC121 

C.CC153 

0.C0183 

8 

0.0C013 

C.CCC24 

0.CCC35 

0.00046 

9 

0.3CC01 

C.0CCC3 

C.CCCC6 

C.C001C 

n 

c.occcc 

c.cccco 

0.CCC01 

C.C0CC2 

u 

c.cocco 

C.CCCCO 

C.COCOO 

12 

O.CCCOC 

C.CCOCC 

13 

O.COGCO 
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APPENDIX  3 (Continued) 


m  -1/14  -1/15  -1/16  -1/17 


2 

0.23862 

C. 23938 

0.24CC5 

0.24063 

3 

0.  1C793 

C.1CSC7 

0.11CC7 

C. 11095 

4 

0.04586 

C. 04694 

C.C4789 

C. 04873 

5 

0.0181C 

C. 01691 

C-C1962 

C.C2C26 

6 

C. 00654 

C.CC7C4 

C.CC749 

C. 00790 

7 

O.OC212 

C.CC239 

C.CC263 

C.C0286 

8 

C.3C06C 

C.C0C72 

C.CCC34 

C.C0C95 

9 

COCO  14 

C.CCL19 

C.CCC24 

C.CC028 

1C 

c.ocro2 

C.CCCC4 

0.CCC06 

C.C0C07 

li 

o.occoc 

C.0CCC1 

C.CCC01 

C.C0CC2 

12 

c.occcc 

c.ccccc 

C.CCCCC 

C.OCCCO 

13 

c.occcc 

c.cccco 

C.CCCCC 

C.CCCCO 

14 

o.ooooc 

c.ccccc 

O.CCCQC 

Q.COOCC 

15 

c.cccco 

O.CCCCC 

o.cocco 

16 

C.CCCCC 

c.ccccc 

17 

c.coccc 

-1/18  -1/19  -1/20  -1/21 


2 

0.24115 

C. 24162 

C.242C4 

G. 24242 

3 

0.11173 

C. 11243 

C. 11306 

0.11363 

4 

0.04948 

C.05C16 

0.C5C76 

0.05131 

5 

C.02°83 

C. 02124 

C.C2181 

G.C2223 

c 

0.0C827 

C.OC861 

C.CCS92 

C.C0920 

7 

C.OC3C7 

C.CC327 

C.CC345 

C.C0361 

3 

C. 00105 

C. 00116 

C.C0125 

C.C0134 

9 

0.0CC33 

C.CCC28 

0.CCC42 

C.C0046 

i: 

C.CC0C9 

C.CCC11 

C.CCC13 

C.C0015 

H 

0.CC0C2 

C.CCCC3 

C.CCC04 

C.CCC04 

12 

c.cccc : 

C.0CCC1 

C.CCC01 

0.C0001 

13 

O.CCCCC 

C.COCCC 

C.CCCCC 

C.OOCCO 

14 

C.CCCCC 

C.CCCCC 

C.CCCCC 

C.OOCCO 

15 

C.OCCCC 

C  .CCCCC 

C.CCCCO 

C.COCCO 

16 

O.OCCOC 

C.CCCCO 

C.CCCCC 

O.COOCO 

17 

C.0C03C 

C.CCCCC 

C.CCCCO 

C.OOCCO 

18 

n.OCOCC 

C.CCCCC 

C-CCCOC 

O.OCCCO 

19 

C.CCCCC 

C.CCCCC 

C.COCCO 

2C 

O.COCOO 

c cocoa 

2  I 

c.occco 

i 

2 
I 
4 
5 

6 
7 
d 

9 


i 
2 

— • 

5 
6 

7 
6 

:: 
u 

12 
13 


APPENDIX  4 

P    (Y=r) 
n 


n 


1  0.5CCC  C.1755                      C.0489  0.0115 

2  l.CCCC  0.8245                      0.5CCC  0.2348 

3  l.CCCC  C.9511  0.7652 

4  l.CCCC  C.9385 

5  l.OOCO 


10 


0.C024 

C.CCC4 

C.CCC1 

c.occo 

0.0904 

0.0298 

C.CC87 

0.0023 

0.5CCC 

C.2695 

C.1234 

0.0493 

0.9C96 

C.73C5 

C.5CCC 

0.2934 

0.9976 

3.97C2 

C.8766 

0.7C66 

l.CCCC 

0.9996 

C.9913 

0.95C7 

l.CCCC 

C.9999 

0.9977 

l.CCCC 

l.COCC 
l.OCCO 

11  12  13  14 


O.CCCC 

c.cccc 

C.CCOC 

O.OCCO 

C.C0C5 

c.ccci 

C.CCCC 

O.OQCC 

0.G175 

0.0C57 

C.CC17 

0.0005 

0.  1499 

C.C678 

C.0277 

0.0103 

0.5C0C 

C.31U 

C.1717 

0.0851 

C.8  5C1 

0.6669 

C.5CCC 

0.3245 

C.98  2  5 

0.9322 

C.8283 

0.6752 

0.9995 

0.9943 

C.9723 

0.9149 

l.CCCC 

C.9999 

C.9983 

0.9397 

l.CCCC 

l.CCCC 

l.CCCC 

0.9995 

l.CCCC 

l.CCOC 

l.OOCO 

l.CCOC 

l.OOCO 
l.COCC 

96 
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APPENDIX  4    (Continued) 


i 
2 
3 
4 

5 
6 

a 

9 
10 

n 
12 

13 
14 
15 
16 
17 


15  16  17  18 


o.occc 

O.CCCO 

O.CCCC 

O.OCCO 

o.ccoc 

C.CCCC 

C.CCCC 

O.OCCO 

O.C0C1 

c.cccc 

C.CCOC 

O.CCCO 

0.0035 

c.ccu 

C.CC03 

0.0001 

0.0  38  2 

0.0158 

0.CC6C 

0.0022 

0.1899 

0.1CC8 

C.C49C 

C.0P2C 

f\      c  **  r  r 

• 

C.3259 

C.2C55 

0.1152 

0.8101 

0.6641 

C.5CCC 

0.3452 

0.9fcl7 

0.8992 

C.7945 

0.6548 

0.9965 

C.9E42 

G.951C 

0.8848 

0.9  99  9 

0.9989 

0 . 9  94  C 

0.9780 

1.C00C 

l.CCCO 

G.9997 

0.9978 

l.CCOC 

1* vlwC 

l.CCCC 

0.9999 

1.CC0C 

l.OCCC 

l.CCOC 

l.OCOO 

l.CCCC 

l.CCCC 

l.OCCC 

l.CCOC 

l.OCOO 
i.OOCC 

19  2C  21  22 

1  O.CCCC  C.OCCC  C.CCOC  O.OCCO 

2  O.CCOC  C.CCCC  C.CCOC  O.OCCO 

3  O.CCCC  C.CCCC  C.CCCC  O.OCCO 

4  O.CCCC  O.CCCO  C.CCCC  O.OCCO 

5  0.CCC7  C.CCC2  C.0C01  O.OCOi 

6  0.CC92  C.0C26  C.CC13  0.0005 

7  0.0595  0.C266  C.C128  0.0C54 

8  0.219C  C.1282  C.C697  C.0353 

9  0.5C0C  C.353C  C.2308  0.1401 

10  0.781C  C.647C  C.5C0C  0.3598 

11  0.94C5  C.8718  C.7692  0.6402 

12  0.99C6  C.9714  C.93C3  0.8599 

13  0.9992  0.9964  C.9872  0.9647 
i*  l.OOOC  0.9998  C.9987  0.9946 

15  l.COOC  l.CCCC  0.9989  0.9995 

16  l.OCCC  l.CCCO  l.CCCC  0.9999 

17  l.CCCC  l.CCCO  l.CCCO  l.OCOO 
13  l.CCOC  l.CCCO  l.CCOC  l.CCOO 
I9  l.CCCC  l.CCCC  l.OCOO 
2C  l.CCCC  l.OCCO 
21  1.0000 
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